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ABSTRACT 

Chaff is the least expensive and simplest soft kill measure 
against anti-ship missiles. In this thesis we study and simulate 
the deployment of chaff in centroid mode for anti-ship missile 
defense. The idea is to fire chaff of radar cross section greater 
than that of the ship, cause a shift of centroid and achieve break 
lock by moving the ship out of the radar resolution cell of the 
missile. After firing the chaff in an optimum direction, the ship 
is made to maneuver in two steps; first to make an earliest 
possible break lock and then to maximize its displacement from the 
missile line of flight. Since the interaction scenario is time 
variant, modelling has been structured on a relative coordinate 
system. The algorithms to implement these objectives are based on 
values extracted from this coordinate system centered on the chaff, 
mobile with the wind and oriented towards the missile. 

The mathematical model is translated into a simulation 
program. We consider it to be a good tool for insight and 
understanding and feel that it has a promise of application both 
in design work and at sea in the passive point defense problem. 
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THESIS DISCLAIMER 



The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of computational 
and logic errors, they cannot be considered validated. Any 
application of these programs without additional verification is 
at the risk of the user. 
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I . INTRODUCTION 



A. IMPORTANCE OF DEFENSE AGAINST ANTI-SHIP MISSILES 

The maritime combat scenario of today is filled with complex 
multi-dimensional threats for a man of war. Conflicts in the past 
quarter century or so have shown that warships have grown more 
vulnerable to aerial threats, in particular to anti-ship missiles. 
We have the following prominent episodes of recent times to endorse 
and highlight this fact: 

• An Israeli destroyer EILAT was sunk near Port Said in the 1967 
Arab-Israeli war by Egypt with Styx missile fired by Soviet built 
Komar class boats. 

• The Pakistan Navy destroyer KHAIBAR and minesweeper MUHAFIZ 
were sunk and logistic ship DACCA was damaged by salvoes of Styx 
missiles fired by Ossa class missile boats, which are also Soviet 
built, during the 1971 Indo-Pakistan war. 

• Gabriel missiles launched by the Israeli Navy inflicted heavy 
damage to the Syrian Navy resulting in sinking of a number of 
missile boats during the 1973 Arab-Israeli war. 

• The Royal Navy suffered a humiliating blow when a number of 
its front line ships were sunk and damaged by shore and air 
launched Exocet missiles fired by the Argentinian Navy during the 
Falklands war in 1982 [Ref. 10]. 

• US ship and air launched Harpoon missiles destroyed a number 
of Libyan missile boats and other craft in March 1986 during a 
conflict with Libya in the Gulf of Sidra. 

• In April 1988 again. Harpoon and Standard missiles fired by 
the US Navy caused heavy damage to Iranian Navy vessels in the 
Persian Gulf. 



B. APPROACHES TO SHIP DEFENSE 

All major and minor maritime forces of the world recognize the 
vulnerability of surface ships to this serious threat today. 
However, since the surface fleets are foreseen to be the mainstay 
of most navies for a long time to come, their defense against anti- 
ship missiles has received tremendous emphasis. Larger navies who 
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have access to sophisticated technological and vast monetary 
resources tie their belief to the methodology of layered defense 
or defense in depth. Such a defense methodology would probably 
include offensive long-range aerial strikes against enemy storage 
depots and platforms, combat air patrols, standoff missiles and a 
variety of other active or passive tactics. Even these elaborate 
measures fail to provide an absolute defense. Smaller navies who 
do not have such luxuries available at their disposal must live 
with lesser and cheaper methods. They have to rely on early 
detection of the enemy and in most cases early detection of the 
missile. Such a restriction imposes the use of short-range active 
or passive measures. These measures include active ECM, 
conventional guns, point defense systems or close in weapon systems 
(CIWS) such as Vulcan Phalanx/Breda and chaff. Ideally one could 
achieve most effective results by having all of them working in 
unison in an intelligent manner. However, availability of active 
ECM and CIWS involves sizable finances and in many cases their 
availability is dependent upon politics. Chaff, on the other hand, 
presents itself to be a cheap and simple alternative. 

C. MODES OF CHAFF DEPLOYMENT 

Chaff has been used effectively in a number of military 
conflicts since the second World War. Its flexibility enables it 
to be used at various stages of conflict in many modes. The 
discussion here will be restricted to a confrontation between an 
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anti-ship missile and a surface unit at sea. Chaff, in such a 
situation can be used in the following modes [Ref. 13 :p. 338]: 

1. Confusion 

In confusion mode, chaff is fired at long range (a few 
kms) when the ship is outside the range of enemy search and target 
indication radars. This tactic provides a large number of false 
targets that are indistinguishable from the real ship. 

2 . Distraction/Dilution 

When the ship has been detected by the enemy and a 
missile attack is imminent, chaff is fired at shorter ranges (1-2 
kms) in order to mislead the missile radar when it starts 
transmitting. The chaff cloud's radar cross section need only be 
comparable to the ship in this mode. 

3 . Dump 

This mode is also called range gate pull off/stealing. 
Once the missile is locked on, chaff is fired at short ranges (less 
than 1 km) . In this mode the ship's ECM equipment is also used in 
conjunction with the chaff. The ship's jammer misleads the missile 
on to the chaff and when jamming ceases, the missile homes onto the 
chaff decoy. In other words, break lock is achieved with the help 
of the jammer. 

4 . Centroid 

In this mode, chaff is deployed at very close ranges. 
The exact figure depends on the philosophy of actual deployment. 
Some experts suggest firing chaff at a range of 25 to 50 m so that 
its echo merges with that of the ship and the missile homes onto 
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the missile is lured away avoiding a direct hit thus minimizing 
damage. This is a last ditch method as break lock is not the 
primary goal. 

D. RECOMMENDED SOLDTION 

Although reaction times are shrinking with sophisticated 
missiles entering the battle arena, intelligent deployment of chaff 
has been known to present an effective counter. Ideally, chaff 
could be most effective if integrated with other forms of active 
and passive anti-ship missile defenses. However, if monetary 
constraints preclude such a luxury, we consider that chaff provides 
a good alternative to thwart most missile threats of today. It can 
be argued that recent advances in technology have enabled 
production of smart and sophisticated missile radars which have the 
ability to reject chaff clutter by sensing the difference of its 
parameters from those of the ship, thereby rendering it 
ineffective. Nevertheless, there is always another counter to a 
counter. Chaff systems are also getting smarter with added 
features. These features include chaff's maneuvering capability 
with the help of a self propulsion system inside the chaff and 
expendable ECM devices with integrated processing facility. Active 
decoy developments such as Nulka and Winnin are examples of such 
chaffs. 

E. AIM OF THE THESIS AND RELATED WORK 

The purpose of this work is to simulate centroid mode of chaff 
deployment, fire the chaff at close range (100 to 400 m) , achieve 
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break lock and maximize miss distance by intelligent ship maneuvers 
assuming a supporting jammer is not available. Work has also been 
done in this field by David Nitzan and Chris Nygard [Refs. 22 & 23] 
The mathematical modelling and simulation done in this thesis has 
been developed as an extension to their work. 
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II. KINEMATICS OF THE PROBLEM 



A. MISSILE DETECTION 

One of the most formidable and decisive problems in the field 
of anti-missile defense is early detection of the missile. Most of 
today's anti-ship missiles travel at speeds close to MACH 1 and fly 
very low to the surface of the sea. This sea skimming capability 
coupled with the state-of-the-art aerodynamic characteristics, 
makes the task of missile detection even more difficult for the 
Operations Room of the ship. For vessels operating without any air 
support, reliance has to be placed on sensors available on board. 
The two key sensors of help in this situation are radar and ESM 
equipment. ESM can be coupled with radar to confirm detection made 
difficult by background clutter at low altitudes. However, modern 
missiles commence radiating for their terminal homing phase quite 
late in their flight, thus limiting the ESM capability to ranges 
less than the radiation commencement zone. 

The primary sensor, as a matter of principle, should be radar. 
Nonetheless, because of low flying profiles of the missile and non- 
availability of state-of-the-art radars to smaller navies, ESM 
takes over the primary sensory role. This is especially true for 
the sea skimmers. The range at which the missile commences its 
terminal phase has a very large variability, depending on the type 
and generation of the missile technology. Based on available 
unclassified data, this range varies from around 3 miles to 20 
miles. With an average speed of 600 knots, at which most missiles 
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travel, we get a reaction time ranging from 18 seconds to about 2 
minutes. It is quite obvious that reaction time is the key factor 
in today's warfare. With the missiles becoming supersonic and their 
terminal phases rapidly shortening, reaction time is increasingly 
diminishing. Table 1 shows reaction times for different initial 
detection ranges and is a useful indicator to emphasize this point. 
Assuming the missile speed to be 600 knots, the following data has 
been tabulated: 

TABLE 1 

MISSILE RANGE AND REACTION TIME 



S.NO 


Missile Ranaefvds) 


Reaction Time (secs) 


1. 


20,000 


60 


2 . 


15,000 


45 


3. 


14,000 


42 


4 . 


13 , 000 


39 


5. 


12,000 


36 


6. 


11,000 


33 


7. 


10,000 


30 


8. 


9000 


27 


9. 


8000 


24 


10. 


5000 


15 


11. 


3000 


9 


12 . 


1000 


TOO LATE 



B. CHAFF DEPLOYMENT IN CENTROID MODE 

There are two philosophies associated with centroid or 
seduction mode tactic; the first is to fire the chaff at very close 
ranges, i.e., 25 to 50 meters [Ref. 5], within the radar resolution 
cell of the missile. The missile takes the two echoes to be one and 
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locks on to the centroid of the two. The ship then moves away 
briskly from the chaff, depending on the wind direction leaving the 
missile to home on the chaff. Depending on all the variables 
involved, the ship can either break lock if it is out of the 
resolution cell in time or else can maximize displacement from the 
missile impact point. In this way, the ship can either escape 
damage altogether or at least, minimize it. Real time data from 
relevant conflicts is not available to seek a recommended figure 
for safe distance of the ship from the point of impact. However, 
at least 50 to 100 meters is considered to be a safe distance for 
minimal damage to the ship. 

The second philosophy behind utilization of chaff in centroid 
mode is to fire the chaff within the resolution cell on a bearing 
related to the wind direction and speed. Chaff is placed as far as 
possible from the ship towards the boundaries of the resolution 
cell [Ref. 13], depending on the initial detection range of the 
missile, beam width of the missile radar and the wind parameters. 
Once the chaff has bloomed, the missile amends its homing on 
direction to the centroid. In the meantime, concurrent to chaff 
launch, the ship maneuvers its way out of the resolution cell to 
achieve a break lock and maximize distance from the missile impact 
point. Analysis of this philosophy is the aim of this thesis. 

C. RADAR RESOLUTION CELL 

The resolution cell is the central factor of the anti-ship 
missile problem. Basically, the radar has a certain range and 
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bearing resolution based on its effective pulse width and beam 
width [Ref. 3:p. 332]. The area outside these bounds is the limit 
of the radar to distinguish one target from another. If there are 
two targets and the radar is locked onto one of them, the second 
should be outside the boundary to be discriminated as a separate 
target. If the second target also lies within this area, the radar 
will only see them as one target and align itself to their 
centroid, both radially and azimuthally. Figure 1 shows what a 
radar resolution cell looks like. 



Range Gate 




Figure 1. Radar Resolution Cell [Ref 23] 

The limits in both angle and range are defined by the half 
power or 3 dB values. Most missile seekers use this geometry for 
processing purposes. From the missile's point of view, only the RF 
energy inside this cell is considered to be the useful backscatter. 
Any return from the target outside this area is ignored to ensure 
sharp resolution. The cell actually assumes an annular shape, but 
at ranges we are interested in, it can be approximated to be a 
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rectangle without incurring any sizable inaccuracies as shown in 
Figure 2: 



Beamwidth 




APPROXIMATION 



Figure 2. Approximation to the Resolution Cell [Ref. 23] 

The bounds of the cell are also called gates; the one 
perpendicular to missile direction is the angle gate and the one 
parallel to it is the range gate as shown in Figure 3. 

MISSILE 

J 

J RANGE GATE 

◄ ► 

ANGLE GATE 

Figure 3. Range and Angle Gate 
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The center of the cell is the centroid for the missile radar 



to lock on. As the bearing of the missile changes due to movement 
of the ship or centroid, the cell rotates to reorient itself in 
order to maintain orthogonality with the missile line of flight. 
Similarly, due to variation of target range from the missile as it 
closes in, the dimensions of the cell vary in the azimuthal plane. 
The radar resolution cell, therefore, forms an orthogonal 
coordinate system which is varying with time and constantly 
readjusting its bounds. This motion, as mentioned, depends on the 
movement of the centroid and the longitudinal axis of the missile. 

D. EFFECT OF CHAFF ON RADAR RESOLUTION CELL 

The whole idea of using chaff in a favorable manner is to 
shift the centroid of the resolution cell away from the ship, walk 
it off to the chaff and move out of the cell to break lock thereby 
seducing the missile onto the chaff cloud. The chaff, of course, 
has to be fired to bloom inside the resolution cell. If the radar 
cross section of the chaff is comparable to that of the ship, the 
missile is expected to shift its lock to the centroid of the two 
radar cross sections which is approximately the mid-point on the 
line joining the ship and the chaff [Ref. 13:p. 347]. 

Usually, most captains would like to fire chaff with a radar 
cross section considerably larger than that presented by the ship. 
The shift of the centroid would then be more favorable to the ship 
and closer to the chaff cloud. In any case, as soon as the chaff 
blooms, it starts moving with the wind. The aim is to fire the 
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chaff on a bearing that optimally assists the ship for the wind 
vector existing during interaction. The job of the ship then is to 
move out of the resolution cell, in a direction that will take it 
away from the chaff. 

The objective of moving the ship out of the cell is to reduce 
and then eliminate its radar cross section from the bounds of the 
cell. As the ship starts moving away from the chaff, its distance 
from the centroid also increases. As it crosses the cell boundary, 
either through its range or angle gate, its radar cross section 
inside the cell decreases. The effect of this development is two- 
fold: first the ship has moved away from the centroid; second, due 
to its decreased radar cross section inside the cell, the centroid 
moves further towards the chaff. This effect multiplies as more of 
the ship exits the cell and hastens the shift of centroid to the 
chaff. As the ship leaves the cell completely, break lock is 
achieved. The missile now locks onto the chaff and the ship 
endeavours to maximize its distance from the new centroid, i.e., 
the chaff. 

The kinematics are, however, not as simple as they appear to 
be. The most crucial factor at this stage is the wind speed and 
direction in relation to both the missile line of flight and the 
ship's initial course. The wind could be a helping hand or prove 
unfriendly depending on what it is doing during the operation. If 
it is strong and moving parallel to the missile line of flight, 
i.e., perpendicular to angle gate, the chaff will spend very little 
time in the cell. This time might not be enough to enable the ship 
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to exit the cell before the chaff does so. As a result, the missile 
locks on to the ship again. 

As an exception, the ship and the chaff might decide to leave 
the boundaries of the cell simultaneously. In such a situation, it 
is hard to tell what the missile will lock on to. There appears to 
be an equal probability of the missile either locking on to the 
ship or the decoy. However, it first depends on the decision-making 
algorithm of the processor of the missile's radar. It would also 
depend on the relative location of the chaff and the ship with 
respect to the missile. If the processor unit conducts a far-to- 
near processing pattern, the nearer of the two will be selected. 
Of course, the reverse would be true for a near-to-far pattern. 

In summary, the presence of chaff redefines the boundaries of 
the radar resolution cell, reorients it, displaces it and aids the 
ship to get out of troubled waters. The crucial factor is a 
favorable wind vector at the scene of the battle. 
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III. MATHEMATICAL MODELLING 



In order to formulate a simulation program, it is necessary 
to construct a mathematical model of the situation in which the 
chaff is to be deployed. There are a number of parameters that are 
made available to us by the ship, missile, wind and the chaff. Some 
of these parameters are constant and others are variable. The 
variables are classified into two categories; those that are 
presented to us by a particular situation or the input variables, 
and those that we have the ability to control and manipulate. It 
is the latter type of variables that make it possible for us to 
successfully use the chaff. A list of input and control variables 
are given in Tables 2 and 3. The abbreviations for all the 
variables in this chapter are chosen to correspond to those used 
in the simulation program. 

A. SHIP'S EQUATIONS 

As summarized in the tables, the ship's parameters consist of 
its length, breadth, course, speed, acceleration, turn rate and 
radar cross section. 

1. Initial Conditions 

The ship is assumed to be at the center of the reference 
coordinate system for the initial set up at the time the missile 
is detected. This is actually the screen center. The initial 
velocity vector of the ship is determined by its initial course and 
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speed. The following equations describe the ship at the starting 
point: 



XS(h) 


= XCENT 


(3.1) 


YS(h) 


= YCENT 


(3.2) 


XSD(h) 


= VS(h)*Sin{ (CSE(h)*)r/180} 


(3.3) 


YSD(h) 


= VS(h)*Cos{ (CSE(h)*)r/180) 


(3.4) 



TABLE 2 

LIST OF INPUT VARIABLES 



S.No. 


VARIABLE 


ABBREVIATION 


MISSILE 


1. 


Initial Bearing 


BRGM(l) 


2. 


Initial Range 


RANGE 


3 . 


Missile Speed 


VMM 


4 . 


Radar Beamwidth 


e 


5. 


Range Gate 


RGATE 


6. 


Missile Acceleration 


AM 


7. 


Proportional Navigation Constant 


N 


SHIP 


8. 


Course 


CSE(l) 


9. 


Speed 


VS 


10. 


Length 


SL 


11. 


Breadth 


SB 


12 . 


Turn Rate 


TR 


13 . 


Acceleration (Velocity Rate) 


VR 


14 . 


Max Radar Cross Section 


RCSSMAX 


15. 


Min Radar Cross Section 


RCSSMIN 


CHAFF 


16. 


Radius 


CRAD 


17. 


Radar Cross Section 


RCSC 


WIND 


18. 


Wind Direction 


GW 


19. 


Wind Speed 


VW 


MISCELLANEOUS 




20. 


Screen Scale 


SC 



15 



TABLE 3 



CONTROL VARIABLES 



S.No. 


VARIABLE 


ABBREVIATION 


SHIP 


1. 


First Leg of New Course 


CSE2 


2. 


Second Leg of New Course 


CSE3 


3. 


Max Speed 


VSMAX 


CHAFF 


4. 


Bearing to launch chaff 


BRGC 



Where 

XS(h) and YS(h) are coordinates of ship displacement 
XSD(h) and YSD(h) are coordinates of ship velocity 
XCENT and YCENT are coordinates of screen center 
CSE(h) is the ship's course 
2. Ship's Movement 

Having detected the missile, the ship after firing the 
chaff, accelerates and maneuvers out of the radar resolution cell. 
In order to update its position, velocity and acceleration 
components, the following equations are used; 

XS(h+l) = XS(h) + XSD(h) + 0.5 XSDD(h)*h^ (3.5) 

YS(h+l) = YS(h) + YSD(h) + 0.5 YSDD(h)*h^ (3.6) 

where 

XSD(h) = VS*Sin{CSE(h) *7 t/180} (3.7) 

YSD(h) = VS*Cos{CSE(h) *jr/180) (3.8) 

The velocity components of the ship are shown in Figure 4. 
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Figure 4. Ship's Velocity , Components 

Also 

VS(h+l) = VS(h) + VR*h (3.9) 

XSDD(h) = -VS(h) {TR*jr/180) * U * Sin{ CSE (h) *jr/180 ) 

+ VR*Sin{CSE(h) *jr/180) (3.10) 

YSDD(h) = -VS(h) {TR*jr/180} * U * Cos { CSE (h) *7 t/180 ) 

+ VR(h) *Cos{TR*7t/180) (3.11) 

We have supposed that VS can attain a maximum value of 25 
)cnots. The terms used in the above equations are defined as 
follows : 

XSDD(h) and YSDD(h) = Instantaneous coordinates of ship's 
acceleration 

XSD(h) and YSD(h) = Instantaneous coordinates of ship's 
velocity 

VS = Velocity of ship 

VR = Velocity rate or acceleration of the ship 
U = Turn direction (discussed later) 
h = Time step increment 
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Updating the bearing of the missile from the ship is carried 



out using the following equations: 



BRGM(h) = 180/ jr * tan 





(3.12) 



where XM(h+l) and YM(h+l) are updated coordinates of the missile 
position. 

3. Radar Cross Section 

The radar cross section of the ship is a variable 
dependent upon the aspect presented by the ship to the missile 
radar. It is assumed that the missile presents maximum radar cross 
section when viewed from right ahead, right astern and the beams. 
As the aspect changes, the radar cross section decreases to a 
minimum at the bows and the quarters. In the absence of real time 
data, we consider this assumption to be a reasonable approximation. 
It is modelled with the help of the following equation: 

RCSS = RCSSMIN+ | RCSSMAX-RCSSMIN*Cos [ 2 { CSE ( 1) -BRGM(l) )*A/180] | 

(3.13) 

where 

RCSS = Radar cross section of the ship 

RCSSMAX = Maximum radar cross section of the ship 

RCSSMIN = Minimum radar cross section of the ship 

CSE(l) = Initial course of the ship 

BRGM(l) = Initial bearing of the missile 
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Pictorially, the equation produces a graph of the ship's radar 
cross section vs the aspect angle as shown in Figure 5. 




Figure 5. Radar Cross Section Model 



090 
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4. The Ship's Model 

The ship's geometry is approximated by a polygon centered 
on (XS,YS) using its length and breadth. The orientation of the 

ship, pivoted on its center is shown with the help of its course 

with respect to true North as shown in the Figure 6. 

The equations to find its vertices are given as under: 
XSHIP(l) = XS(h) + SL/2 * Sin {CSE(h)) (3.14) 

YSHIP(l) = YS(h) + SL/2 * Cos (CSE(h)} (3.15) 

SHIPS MODEL 




Figure 6. The Ship's Model [Ref. 23] 
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XSHIP(2) = XS(h) 
YSHIP(2) = YS(h) 
XSHIP(3) = XS(h) 
YSHIP(3) = YS(h) 
XSHIP(4) = XS(h) 
YSHIP(4) = YS(h) 
XSHIP(5) = XS(h) 
YSHIP(5) = YS(h) 



+ SL/4 * Sin 



+ SL/4 * Cos 



- SL/2 * Sin 



- SL/2 * Cos 



- SL/2 * Sin 



- SL/2 * Cos 



+ SL/4 * Sin 



+ SL/4 * Cos 



{CSE(h) ) + SB/2 
{CSE(h) } - SB/2 
{CSE(h)} + SB/2 
{CSE(h) } - SB/2 
{CSE(h) } - SB/2 
{CSE(h) ) + SB/2 
{CSE(h) } - SB/2 
{CSE(h) } + SB/2 



* Cos(CSE(h) ) 

(3.16) 

* Sin {CSE(h) } 

(3.17) 

* Cos {CSE(h)} 

(3.18) 

* Sin {CSE(h) } 

(3.19) 

* Cos {CSE(h)} 

(3.20) 

* Sin {CSE(h) } 

(3.21) 

* Cos (CSE(h) ) 

(3.22) 

* Sin {CSE(h) } 

(3.23) 



B. CHAFF EQUATIONS 

The following parameters are related to the chaff: 

• Launch bearing 

• Distance from the ship 

• Radius of the chaff 

• Radar cross section of the chaff 

• Bloom time 

• Velocity vector of chaff after bloom 

1. Initial Conditions 

Since the wind is the fundamental issue here, it is 
appropriate to see its relationship with the chaff. Wind direction 
is, by convention, taken to be the direction it is blowing from. 



21 



the direction in which the chaff moves is its reciprocal. This is 
shown in Figure 7. The unsealed velocity of the wind is VWl and it 
is also the velocity of the chaff. 



Velocity Vector 




Figure 7. The Wind Velocity Vector 

where 

VWl = Wind speed (Unsealed) 

GW = Wind direction 

Velocity vector of the chaff is given by the following 
equation: 

VC = VW 

GC = GW + 180 (3.24) 

Let us now consider the initial positioning of the chaff. In 
order to offset the centroid as far away as possible, we fire the 
chaff so that its center lies on the boundary of the angle gate. 
The angle gate is, however, shrinking as the missile closes in. 
The initial placement of the chaff, therefore, is brought closer 
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to the cell center laterally by one chaff radius in order to make 
it bloom inside the future dimensions of the cell. The calculation 
of this distance is done using the following equation: 

CDIST = R(l) * TAN (THETA * PI/180) - CRAD (3.25) 
The initial position of the chaff round is shown in Figure 8. 



Launching Vector 



North 



Ship 



(XC, YC) 

BRGC V 

CDISV 

CDIST Cos BRGC 




CDIST Sin BRGC 
(XCENT, YCENT) = (XS, YS) 



Figure 8. Initial Position of Chaff [Ref. 23] 

The parameters in the figure are as follows: 

BRGC = Chaff launch bearing 

CDIST = Distance of chaff from the ship 

We can see here that the coordinates of the chaff are: 

XC(h) = XCENT + CDIST * Sin(BRGC * jr/180) (3.26) 

YC(h) = YCENT + CDIST * Cos(BRGC * k / 180 ) (3.27) 

2. Chaff Movement 

The ship is initially in the screen center. As the chaff 
starts blowing with the wind, its coordinates are updated by the 
following equations: 
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XC(h+l) = XC(h) + VC * Sin (GC * ;r/180) 
YC(h+l) = YC(h) + VC * COS (GC * jr/180) 



(3.28) 

(3.29) 



3. Bloom Time Delay 

In order to Maximize the factor of realism, we have taken 
into account the time delay between the firing of the chaff and its 
full bloom. For simplicity, it has been assumed that the time taken 
to launch the chaff is included in the bloom time delay. Although 
various types of chaff have different time delays in this context, 
we consider a bloom time delay of 5 seconds as an appropriate 
average figure [Ref. 13:p. 342]. It also includes the time it takes 
to launch the chaff at its intended initial position. The 
implmentation of this concept is achieved by keeping the chaff on 
the ship for the first 5 iterations. On the sixth second, the chaff 
cloud comes into existence at the desired location. Thereafter it 
is at the mercy of the wind vector and starts drifting. Since the 
ship and the resolution cell are mobile, there is a possibility 
that the chaff may bloom outside the cell as a result of bloom time 
delay. To avoid such an eventuality, the first position of the 
bloomed chaff is offset in accordance with the ship^s movement. 

4. The Chaff Model 

The chaff cloud is approximated by a circle but due to 
convenience of mathematical calculations, it is drawn as a hexagon 
using graphics. The hexagon is centered on the chaff center (XC,YC) 
and is shown in Figure 9. Coordinates for the corners are worked 
out with the help of the following equations: 
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XCF(l) = XC(h) + CRAD 
YCF(l) = YC(h) 



(3.30) 

(3.31) 



XCF(3), YCF(3) XCF(2), YCF(2) 




XCF(2) 


= 


XC(h) 


+ 


CRAD 


★ 


Cos(60' ) 


(3 . 32) 


YCF(2) 


= 


YC(h) 


+ 


CRAD 


★ 


Sin(60“ ) 


(3.33) 


XCF(3) 


= 


XC(h) 


- 


CRAD 


* 


Cos(60° ) 


(3.34) 


YCF(3) 


= 


YC(h) 


•f 


CRAD 


★ 


Sin(60“ ) 


(3 . 35) 


XCF(4) 


= 


XC(h) 


- 


CRAD 






(3 . 36) 


YCF(4) 


= 


YC(h) 










(3.37) 


XCF(5) 


= 


XC(h) 


- 


CRAD 


* 


Cos(60°) 


(3 . 38) 


YCF(5) 


= 


YC(h) 


- 


CRAD 


* 


Sin(60“ ) 


(3.39) 


XCF(6) 


= 


XC(h) 


+ 


CRAD 


* 


Cos(60') 


(3.40) 


YCF(6) 




YC(h) 


- 


CRAD 


★ 


Sin(60' ) 


(3.41) 



C. EQUIVALENT TARGET EQUATIONS 

Before the chaff is fired, the missile radar is locks on to 
the ship and is tracking it. When the chaff is fired and it has 
bloomed inside the radar resolution cell of the missile, the radar 
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senses a new situation. The weighted center of the two, i.e., the 
ship and the chaff, is now seen by the missile as the centroid of 
the target and it starts tracking it. This centroid, also called 
the equivalent target, is determined by the magnitude of the radar 
cross sections of the chaff and the ship as seen by the missile. 
It is calculated by the following equations: 

XET(h.l) = XS(hl «RC||jhj . RCSC(h) 

YET(h.l, = 1811 ii!i : imt,; (3.43) 

The velocity with which this equivalent target moves is calculated 
using the following equations: 

VET = [{XET(h) - XET(h-l)}^ + {YET(h) - YET (h-1) (3.44) 



GET = tan 



YETfh) - YETfh-1) 
XET(h) - XET(h-l) ' 



(3.45) 



This is shown in Figure 10. 



D. MISSILE EQUATIONS 

The missile parameters consist of the following: 

• Range and bearing from the ship 

• Range resolution 

• Missile speed 

• Half power beamwidth of the antenna 

• Proportional navigation constant 

• Acceleration of the missile 

• Turning rate 
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(XET(h-1)},YET(h-1)) 

Figure 10. Equivalent Target Velocity 
1. Initial Conditions 

The first parameters related to the missile that are 
observed by the ship are its initial detection range and bearing 
which are shown in Figure 11. The related equations are as follows: 
XM(i) = XCENT + R(i)* Sin{BRGM(i) * 5 t/180} (3.46) 

YM(i) = YCENT + R(i)* Cos{BRGM(i) * ?r/180) (3.47) 

The direction flight of the missile is the reciprocal of its 

bearing from the ship i.e., 

GM = BRGM(i) + 180 (3.48) 

The initial velocity components of the missile in the X and Y 
directions are shown in Figure 12. The equations for these 

components are as follows: 

XMD(h) = VM * Sin(GM * t/180) 

YMD(h) = VM * Cos(GM * jr/180) 
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(3.49) 

(3.50) 



North 



(XCENT, YCENT) 




Figure 11. Missile Range and Bearing 



North 




Figure 12 . Missile Velocity Components 
2 . Missile Movement 

The coordinates of the missile position as it moves are 
updated as follows: 
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XM(h+l) = XM(h) + XMD(h) * h + 0.5 * XMDD(h) * (3.51) 

YM(h+l) = YM(h) + YMD(h) * h + 0.5 * YMDD(h) * (3.52) 

XMD and YMD are extracted from the equations in the previous 
section whereas, XMDD and YMDD, which represent the components of 
the missile acceleration are shown in Figure 13 and given by the 
the following equations: 

XMDD = AM * Sin(GM) (3.53) 

YMDD = AM * COS(GM) (3.54) 



North 



(XCENT, YCENT) 




(XM, YM) 



Figure 13: Missile Acceleration Components 
The missile direction angle is updated using the following 
equation: 



GM 



tan 



r 




- XM(h 


) 




[HTTy 


- YM (h 


) 



(3.55) 
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This is shown in Figure 14. The X and Y components of the missile 
position vector are updated iteratively. 




XM(h+1) - XM(h) 



Figure 14. Updating Missile Direction 
3 . Target Tracking by the Missile 

In order to maintain track of the target the missile uses 
the principle of proportional navigation constant i.e it maintains 
a constant line of sight with respect to a reference. This angle 
is shown in Figure 15, 
where 



a 



tan 



1 


YET(h 


) - YM(h) 




XET(h 


) - XM(h) 



(3.56) 



In order to maintain a constant line of sight angle, the 



missile guidance system senses the rate of change in o. This is 



given by the following equation: 



d/dt (c7) 



c-rm - VET*Sin( fGET-g) ) - VM*Sin( (GM-p ) ) 

RHrTTil 



(3.57) 
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(YET-YM) 

— REFERENCE 

Figure 15. Target-Missile Line of Sight [Ref. 23] 

where 

RMT(h) = Range of the missile from equivalent target 

In order to cater for this rate of change in the line of 
sight and bring it to zero, the missile needs to turn at a 
corresponding rate. This is called the turn rate or the 
acceleration of the missile and is given by the following equation: 
AM = SIGD * VM * n (3.58) 

where 

n = Missile navigation constant 

E. RELATIVE COORDINATE SYSTEM 
1. Shift of Origin 

In the complex and randomly occurring scenario of 
different ship headings, missile bearings and wind directions, 
movement of the chaff after full bloom can take many paths in terms 
of relative motion. In order to minimize the disadvantageous 
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features of the variables involved, we have decided to have a 
relative coordinate system centered on the chaff round. This is 
accomplished with the help of the following equations: 

XR(h) = XS(h) - XC(h) (3.59) 

YR(h) = YS(h) - YC(h) (3.60) 

where XR(N) and YR(N) are the relative abscissae and ordinates 
respectively. A sketch depicting the shift of the reference system 
is shown in Figure 16. 



i 


L 

„SHIP 

(1,2) 

(-1.1) 


I 

CHAFF (NEW ORIGIN) 




OLD ORIGIN 


(2.1) 

(0.0) 




o . 





Figure 16. Shift of Reference System 
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The original coordinates of the positions of the ship and chaff 
are (1,2) and (2,1) respectively but when we shift the origin to 
the chaff, ship's position using the above mentioned equations 
becomes (-1,1). 

2. Rotation of the System 

The next step is to rotate the new coordinate system to 
the bearing of the missile in the clockwise direction. Such a 
rotation is desirable and handy to calculate the chaff launch 
bearings which are orthogonal to missile direction. This is done 
using the following relationship: 





1 


Cos BRGM(h) 




Sin BRGM(h) | 




XRPfh) 


= 1 






1 


XRfh) 




1 


Sin BRGM(h) 




Cos BRGM(h) 1 


(3.61) 


where XRP 


and XR 


are the new 


and 


previous state 


vectors 


respectively 


. Upon 


multiplication 


with 


the matrix we 


get the 


following equations: 










XRP(h) 


= XR(h) 


* Cos BRGM(h) - 


YR(h) 


* Sin BRGM(h) 


(3.62) 


YRP(h) 


= XR(h) 


* Sin BRGM(h) + 


YR(h) 


* Cos BRGM(h) 


(3.63) 



This mathematical model can be better understood using Figure 17. 
We see that point (0,2) becomes (-1,1.7) when the system is 
rotated. 

3 . Relative Velocity Components 

The next important step is to calculate the relative 
velocity components using the initial values of the ship's velocity 
and the velocity of the wind with which the chaff moves. This is 
shown in Figure 18. 
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Figure 17. Rotated Reference System 
Using Figure 18 we get the following equations: 

XRPD(h) = VS*Sin(CSE(h) -BRGM(h) ) + VW*Sin(GW - BRGM(h)) (3.64) 
YRPD(h) = VS*Cos(CSE(h) -BRGM(N) ) + VW*Cos(GW - BRGM(h)) (3.65) 
where VS and VW are the velocities of ship and wind respectively 
and XRDP(h)/YRDP(h) are the X and Y components of the relative 
velocity respectively. These values of relative velocity components 
are calculated using the initial conditions. 
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MISSILE 




F. PRE-COMPUTATION METHOD 

The random variability of the missile direction, ship's course 
and speed at the time of missile detection along with wind speed 
and direction, present a much more complex problem than is readily 
obvious. It is extremely difficult to obtain the best solution in 
all the possible cases with conventional chaff alone. Efforts are, 
however, made in this work to enhance optimization. In order to 
achieve this, it is considered useful to carry out a rough pre- 
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computation which can give an idea as to what chaff firing bearing 
would entail the largest miss distance. Making use of the relative 
coordinate system, the problem becomes somewhat simpler in 
formulation. As mentioned earlier, the chaff is fired orthogonal 
to the missile bearing, leaving us with two alternatives; one to the 
right or in the positive XRP direction and the other in the 
negative XRP direction i.e., to the left. The distance traveled by 
the chaff in either direction is computed for the entire run of the 
play. A couple of assumptions made in this regard make this a 
rather rough but a pretty good estimate of the situation. The first 
assumption is that the ship's speed is taken to be 25 knots which 
is the maximum speed of the ship assumed in this work rather than 
a gradual increase from the ship's initial speed. The next 
assumption ignores the alteration to the first leg of the new 
courses which is carried out to escape the range gate in minimal 
time. A direct alteration to the second leg from the ship's initial 
course is assumed instead, which is perpendicular to the missile 
direction. The following equations calculate the magnitude of the 
relative motion of the ship with respect to the chaff cloud in 



either direction. 

XRPl = XRPl + XRPD1*T + 0.5*XRPDD*T^ (3.66) 

YRPl = YRPl + YRPD1*T + 0.5*YRPDD*T^ (3.67) 

XRP2 = -XRPl + XRPD1*T + 0.5*XRPDD*T^ (3.68) 

YRP2 = -YRPl + YRPD1*T + 0.5*YRPDD*T^ (3.69) 

where 

XRPDl = XRPDl + XRPDD*T (3.70) 
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YRPDl 



YRPDl + YRPDD*T 



(3.71) 



and 

XRPDD = -ACC*Sin(CSEDIFl) (3.72) 

YRPDD = -ACC*Cos(CSEDIFl) (3.73) 

where 

XRPl/YRPl = Distance traveled by chaff in X/Y direction 
XRPDl/YRPDl = Velocity component of chaff in X/Y direction 
XRPDD/YRPDD = Acceleration component of chaff in X/Y direction 
ACC = Acceleration of the ship 

CSEDIFl = Difference in the ship's initial course and the new 
course 

This is shown in Figure 19. 

G. RANGES AND BEARINGS 
1. Ranges 

There are three different ranges of interest shown in 
Figure 20 that are required for subsequent calculations. The range 
between the missile and the ship center is given by: 

R(h) = [ (XS(h) - XM(h))^ + (YS(h) - YM(h) (3.73) 

Similarly the range between the missile and the chaff center is 
given by: 

RCM(h) = [ (XC(h) - XM(h) + (YC(h) - YM(h) (3.74) 

and the range between the missile and the centroid or the 
equivalent target is: 

RMT(h) = [ (XET(h) - XM(h))^ + (YET(h) - YM(h) (3.75) 
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MISSILE 




Figure 19 . Pre- Computed Distances 
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2 . 



Bearings 

Like the ranges mentioned in the previous section, there 
are also some bearings that are necessary for subsequent 
calculations. These bearings are shown in Figure 21. The bearing 
of the equivalent target from the missile is given by: 



BRGMT(h) = tan*’ 




y 



(3.76) 




Figure 21. Various Bearings from the Missile 



39 



The reciprocal of this bearing, i.e., the bearing of the missile 
from the equivalent target is : 



-1 



BRGTM(h) = tan ' < 



XETfh) - XMrh^ 
VETihj - VM(hj 



(3.77) 



The bearing of the ship from the missile is: 



BRGMS(h) = tan 



-1 



YS(h) - YMfh) 

Xs(h) - XMlhj 



(3.78) 



Lastly, the bearing of the chaff from the missile is: 



BRGMS(h) 



tan 





(3.79) 



H. RANGE AND ANGLE GATES 

The range gate of the missile is approximately a function of 
its pulse width and is an input variable to the problem. The angle 
gate, on the other hand, is dependent upon the distance of the 
missile from the target and the half power beamwidth of the missile 
radar antenna. This is calculated using the following equation: 

AG = RMT(h) * Tan(©B) (3.80) 
We also need to calculate the distances of the ship and the chaff 
from the center of the radar resolution cell, which is the 
equivalent target or the centroid of the two radar cross sections. 
The position of the ship with respect to the radar resolution cell 
center along the range and angle gate directions, is calculated as 
follows ; 
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AGS = R(h) * Sin I BRGMS(h) - BRGMT(h) 1 (3.81) 

RGS = I RMT(h) - R(h) * Cos(BRGMS(h) - BRGMT(h)} | (3.82) 

This is shown in Figure 22. Similarly, the position of the chaff 
is given by the following equations: 

AGC = RMC(h) * Sin I BRGMC(h) - BRGMT(h) | (3.83) 



Resolution 

Cell 




Missile 
(XM, YM) 

Figure 22. Ship's Position Relative to Cell Center 
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RGC = I RMT(h) - RMC(h) * Cos{BRGMC(h) - BRGMT(h)) | 
This is shown in Figure 23. 



(3.84) 



I. RADAR RESOLUTION CELL MODEL 

The radar resolution cell is constructed with the help of the 
range and angle gates of the missile radar. It is centered on the 




Figure 23. Chaff Position Relative to Cell Center 



42 



equivalent or the centroid and is approximated to be a rectangle 
as discussed in Chapter II. The cell is oriented towards the 
missile bearing. As mentioned earlier, the range gate of the 
missile is constant. The angle gate, on the other hand, is 
dependent on the half power beamwidth of the antenna and range of 
the missile so it varies as the missile closes in. The cell is, 
therefore, a time variant rectangle increasing in size as the 
missile's range decreases and reorienting to maintain orthogonality 
with the missile direction. 

The depth of the resolution cell is twice that of the range 
gate and the width is twice that of the angle gate. This is shown 
in Figure 24. 
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The edges along the width of the center line are calculated 



as follows: 

XWl = XET(h) + AG * Cos{ (BRGM(h) ) (3.85) 

YWl = YET(h) - AG * Sin( (BRGM(h) } (3.86) 

XW2 = XET(h) - AG * Cos( (BRGM(h) } (3.87) 

YW2 = YET(h) + AG * Sin( (BRGM(h) } (3.88) 

The corners of the cell are calculated as follows: 

XCELL(l) = XWl + RG * Sin { (BRGM (h) } (3.89) 

YCELL(l) = YWl + RG * Cos{ (BRGM(h) ) (3.90) 

XCELL(2) = XWl - RG * Sin{ (BRGM(h) ) (3.91) 

YCELL(2) = YWl - RG * Cos{ (BRGM(h) ) (3.92) 

XCELL(3) = XW2 - RG * Sin{ (BRGM(h) } (3.93) 

YCELL(3) = YW2 - RG * Cos { (BRGM (h) ) (3.94) 

XCELL(4) = XW2 + RG * Sin{ (BRGM(h) ) (3.95) 

YCELL(l) = YW2 + RG * Cos { (BRGM (h) ) (3.96) 



J. PORTION OF THE SHIP INSIDE THE RESOLUTION CELL 
1. Ship Fully Inside the Cell 

The entire radar cross section of the ship lies inside 
the resolution cell initially. After firing the chaff, the ship 
attempts to maneuver out of the cell. In the process, the center 
of the cell shifts to the equivalent target. The new situation 
alters the radar cross section of the ship as seen by the missile 
in order to calculate the centroid. This cross section depends on 
the portion of the ship inside the cell and is a function of the 
following parameters: 
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• Radar resolution cell center 

• Size of the cell 

• Relative location of the ship center 

• Ship's heading 

The RCS of the ship is assumed to be the maximum presented by 
the ship at a particular aspect angle when it is fully within the 
bounds of the cell as shown in Figure 25. 
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Figure 25. Limit for Maximum RCS of the Ship 
in the Resolution Cell [Ref. 23] 

2. Ship Partially Inside the Cell 

When some part of the ship is outside the cell, we have 
a partial radar cross section case. In order to calculate it, we 
need to shift our reference system from true North to the 
horizontal. This shift is necessitated by the fact that the 
computer calculates angles from the horizontal reference as shown 
in Figure 26. 
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Figure 26. Angle Translation to Fortran 23] 

A conversion is, therefore, undertaken to cater for this new 
reference system. We define a new angle 6, as shown in Figure 27 
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Figure 27. True North to Fortran Conversion [Ref. 23] 
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between the new reference and the direction from which the missile 
is approaching the ship. The aspect angle of the ship is called 
'Saa' and it is calculated using 6. It is defined as the angle 
between the ship's heading and 6 as shown in Figure 28. 



Ship Center 



Resolution 
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Figure 28. Aspect Angle [Ref. 23] 

The ship is approximated as a rectangle for ease of 
calculation. Also shown in this figure are two more parameters; 
namely Sir and S2r. The subscript 'r' merely represents the fact 
that the ship is exiting through the range gate boundary of the 
cell. If it exits via the angle gate boundary, then we replace it 
with subscript 'a'. S2r is the shortest distance between the range 
gate boundary and the ship center and it is given by: 

S2r = I RG - RGS I (3.97) 
Sir represents the portion of half of the ship's length that is 
inside the cell and is calculated with the help of the following 
equation: 
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sir 



Cos (Saa) 



(3.98) 



In the case shown in Figure 28, where the ship center is 
inside the cell, the portion of the ship's RCS in the cell is given 
by: 

PARTS = 4“ + IMa <3.99) 



where 

SI = Sir or Sla depending upon how the ship exits the cell 
Sarea = SL*SB 

In case the ship center lies outside the resolution cell, the 
calculation is as follows: 

PARTS = 4- - 11^ (3.100) 



If the ship is exiting a corner (leaving the range and angle 
gate boundaries simultaneously) , then the amount of ship's RCS in 
the cell is approximated by using the lesser of the two values 
obtained from Sir or Sla. 

K. PORTION OF CHAFF INSIDE THE CELL 
1. Chaff Fully Inside the Cell 

The chaff is fired to bloom inside the radar resolution 
cell. After it blooms the wind takes over and the chaff cloud 
starts moving with the wind vector. Since the cell is centered on 



48 



the equivalent target, some portion of the chaff might, dependent 
on the situation, spill out of the cell. In order to get an 
estimate of the location of the center of the cell, it is important 
to know how much of chaff is inside the cell. This depends on the 
following parameters: 

• Relative location of chaff center 

• Size of chaff 

• Location of resolution cell center 

• Size of the cell 

The chaff lies completely within the cell when the distance 
of its center from the cell bounds is at least equal to the radius 
of the chaff. This is shown in Figure 29. 




Figure 29. Limit of Maximum Chaff RCS in the Cell [Ref. 23] 

2. Chaff Partially Inside the Cell: 

We define the following parameters to assist in 
calculation of the chaff's RCS inside the cell: 

Carea = t * (CRAD)^ (3.101) 
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where 



CRAD = Chaff radius 

Rlr = I RGC - RG I (3.102) 

Rla = I AGC - AG I (3.103) 

If the chaff is located at the range gate boundary, the 
situation will be as shown in Figure 30. In this particular case, 
more than half the chaff lies outside the bounds of the cell. 




Figure 30. Portion of Chaff RCS in the Cell [Ref. 23] 

If we expand the view of interest, as shown in Figure 31, we can 
calculate a few more parameters that would assist us in finding 
chaff RCS inside the cell. 

We see that: 



A2r = Cos 



Rlr 

CRAD 



(3.104) 
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(3 . 105) 



The sector bounded by twice this angle is given by: 
SECTr = 1/2 (2*A2r) CRAD^ 




Figure 31. Partial Chaff RCS in the Cell Blown Up [Ref. 23] 

Now 

R2r = CRAD * Sin(A2r) (3.106) 

The area of the two triangles, which is also the portion of the 
sector outside the radar resolution cell is: 

TRIr = 2 (I/2*Rlr*R2r) = Rlr * R2r (3.107) 

In case the chaff exits the cell through the angle gate, 
the same procedure is applied using subscript Rla. If both the 
gates are cut off simultaneously, then we use both procedures to 
calculate the area or the chaff RCS inside the cell. All three 
possible cases are shown in Figure 32. 

When Rlr or Rla are zero individually, the portion inside 
the cell is one half of the total. However, when they are both zero 
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simultaneously, as in the corner case, then it is only one quarter 
of the total . 



PARTc = 0.5 




Figure 32. Various Cases of Partial Chaff RCS [Ref. 23] 

If RG >RGC, then most of the chaff cloud lies inside the 
resolution cell. The portion of chaff in the cell is calculated as 
follows : 



PARTC = 



CAREA - SECTr + TRIr 
CAREA 



(3.108) 



IF RG < RGC, then the RCS in the cell is less than half and is 
given by: 



PAPT'c ■“ SECTr " TRIr 
“ CAREA 



(3 . 109) 



When the chaff exits the cell from a corner we have four 
situations : 

1 . RGC ^ RG and AGC ^AG 

This is shown in the Figure 33 along with the actual and 
approximated radar cross sections. 
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Figure 33. Chaff Exiting the Cell Comer [Ref. 23] 
2 * RGC K RG and AGC ^ AG 
This situation is shown in Figure 34. 
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Figure 34. RGC.LT.RG & AGC.GT.AG [Ref. 23] 
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3. RGC> RG and AGC < AG. This is shown in Figure 35. 
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Figure 35. RGC.GT.RG & AGC.lt. AG [Ref. 23] 

4 • RGC ^ RG and AGC ^ AG 

In this case the larger of the two areas produced in the cell by 
RGC or AGC is approximated to be the chaff RCS of interest as shown 
in Figure 36. Less than a quarter of the chaff is inside the cell 
and it is taken to be one half the area cut off by the lesser AGC | 
or RGC . 



RCS Approximation 





Figure 36. RGC.GT.RG & AGC.GT.AG [Ref. 23] 
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L. CHAFF LAUNCH BEARING AND SHIP MANEUVERS 

1. General Description 

In order to maximize time to break lock and miss 
distance, a double maneuver by the target ship is simulated. The 
first maneuver is to get out of the radar resolution cell of the 
missile via the shortest route, i.e, either towards the missile or 
directly away from it. The next maneuver is to alter to a course 
perpendicular to the missile direction when break lock is achieved 
and the ship is outside the resolution cell of the missile's radar. 

2. Chaff Launch Bearing 

First of all, the bearing to launch the chaff is worked 
out. It is basically desired to launch the chaff at right angles 
to the missile direction. The two possible true bearings are 
calculated as follows: 

BRGCl = BRGM(l) + 90 (3.110) 
BRGC2 = BRGM(l) - 90 (3.111) 

The decision regarding the chaff firing bearing is 
dependent upon the results of the pre-computation in order to 
optimize the miss distance. If the ship travels more along the 
positive XRP direction, we fire the chaff towards BRGCl and vice 
versa. 

In the particular case shown in Figure 37, BRGC is equal 
to BRGC2. It must be kept in mind that the missile is shown to come 
from the top because of the rotated reference system and could in 
actuality be approaching from any direction in the azimuth. 
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Figure 37. Pre-Computation for Calculating BRGC 
3. First Maneuver 

The next object Is to calculate the first course to steer 
to get out of the radar resolution cell after having fired the 
chaff. This situation Is shown In Figure 38. 

The choice to select one of the course options is based 
on the relative movement of the ship along the vertical axis, l.e., 
YRP direction which is oriented along the missile line of flight. 
The following algorithm in the program enables us to decide what 
course option to take: 
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If YRP(N) > 0 then CSE2 or the first leg is into the 

missile . 

If YRP(N) < 0 then CSE2 is away from the missile. 



I MISSILE 



CSEOPT(I) 




I CSEOPT(2) 



Figure 38. Course to Steer 
4 . Second Maneuver 

Once the ship is out of the radar resolution cell we 
calculate the next course to steer to maximize miss distance. This 
is initiated when one of the following conditions occur: 



a) YRP(N) > RG 

b) XRP(N) > AG 

In either case, the object is to ascertain whether the 
ship is out of the radar resolution cell, implying that break lock 
has occurred. However, the subroutine for calculating the 
parameters for the second maneuver only calculates CSE3 and the 
turn direction. The actual implementation of affecting this course 
alteration based on the two conditions already mentioned, is 
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achieved in the main program. Since the second leg of the new 
course, CSE3, is at right angles to the first leg, CSE2, there are 
again two choices presented as shown in Figure 39. 

It has been supposed here that CSE2 is equal to CSE0PT2, 
referring to Figure 38. 

Moreover, 

CSE0PT(3) = BRGM(l) + 90 (3.112) 

CSEOPT(4) = BRGM(l),- 90 (3.113) 
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Figure 39. Second Leg of the Course 
The selection of CSE3 from these two options is related 
to the bearing in which the chaff is launched and is usually in its 
opposite direction. In the situation shown in Figure 39, 

CSE3 = CSE0PT(4) 

The direction to turn is calculated in a manner similar 
to the one described in Section 3. 
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IV. DESCRIPTION OF PROGRAM 



The simulation program is written in FORTRAN-77 (Microsoft 
version 3.2). In order to show the run of play at various times, 
the program uses a compatible Plot-88 (1984) graphics package. 
A. INITIAL CONDITIONS 

The initial conditions are set up through an interactive 
exchange with the operator in order to input all the important 
variables. The following variables are used as inputs and they have 
been assigned values that are considered reasonably standard but 
they can be altered through a simple modification in the program. 

TABLE 4 

INPUT VARIABLES WITH TYPICAL VALDES 



S.NO 


Variadile 


Description 


Value 


1. 


CRAD 


Radius of Chaff 


50 yards 


2. 


SL 


Ship Length 


200 mtrs 


3 . 


SB 


Ship Breadth 


20 meter 


4 . 


VSMAX 


Max Ship Velocity 


25 knots 


5. 


N 


Missile Nav Constant 


3 


6. 


VS(1) 


Initial Ship Vel 


20 knots 


7. 


VMM 


Missile Velocity 


600 knot 


8. 


TR 


Turn Rate of the Ship 


5‘/s 


9. 


VR 


Ship's Acceleration 


0.5 kts/s 


10. 


AM 


Missile Acceleration 


30 yards 


11. 


CDIST 


Range of Chaff 


200 yards 
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A scaling factor is included in the input variables to enable 
actual ranges, velocities and accelerations to be depicted on the 
screen of the monitor. The screen center (4. 0,3.0) is used as a 
reference with (0.0, 0.0) being the lower left hand comer and 
(8.0, 6.0) the upper right hand comer of the screen. 

B. RUN TIME 

The program is incremented in steps of one second intervals 
until the missile gets to within 500 yards from the ship. 
Thereafter, to show details of action, we use 0.1 second intervals. 
The total time to go is found by the following relationship: 



TTG 



RANGE 

VMM 



(4.1) 



where 

RANGE = Initial Range of the Missile 

The total program step intervals including 0.1 second steps are 
given by: 



ITTG = Integer TTG + ' 



(4.2) 



C. MAIN BODY OF THE PROGRAM 

The main body of the program iterates through an array with 
a maximum of 60 iterations allowed. The values of missile, chaff 
and ship parameters are updated along with other related parameters 
in each time step. The related parameters include gates, angles. 
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equivalent target parameters and radar cross sections, etc. The 
missile radar seeker is turned off when its range from the 
equivalent target is 250 yards resulting in zero missile 
acceleration and from thereon, the missile motion is essentially 
ballistic. When the missile has flown through the chaff cloud, its 
range starts increasing again and the main loop is completed on 
that iteration. If the missile passes within 50 yards of the ship, 
it is considered to have hit the ship. The implementation of 
alteration of courses to maneuver the ship is also achieved in the 
main body of the program using the decision algorithms mentioned 
in Section ' h * . 

D. USER INTERACTION AND OUTPUT 

After having initialized all the variables, the user is 
presented with two options; either to let the program calculate the 
control parameters by resorting to auto mode or to manually inject 
these parameters. The parameters consist of chaff firing bearing, 
both legs of the courses to be steered to maneuver the ship and the 
turn directions. Subsequent to this stage, the user is presented 
with a picture comprising of initial conditions as shown in Figure 
40. The essential parameters are listed on the picture to its 
right and at the bottom. Thereafter, the user is presented with 
another menu of the following options: 

1. Option One 

This option draws the missile, ship and chaff paths and 
gives the miss distance at the bottom as shown in Figure 41. 
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Figure 40. Initial Conditions 
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Figure 41. Paths Traversed by Missile, Ship & Chaff 
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MINIMUM RRNGE TO SHIP 836.1 YARDS 



2 . Option Two 

This option provides a pictorial representation quite 
similar to the initial conditions picture and we see the situation 
in the resolution cell in 5 second increments. It depicts the run 



of play showing successive positions of the missile, chaff and the 



ship. An example is shown in Figure 42. 
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Figure 42. Instantaneous Position of the 
Resolution Cell 
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3. Option Three 

Option three facilitates changes in the input scenario 
parameters with the help of a menu. It is incorporated to provide 
flexibility to the user to make any alterations in the input 
parameters as desired. The parameters can be changed by choice 
either individually or in totality. The menu of options is as 
follows: 



TABLE 5 

CHANGEABLE PARAMETERS 



• MISSILE BEARING 

• MISSILE RANGE 

• SHIP'S HEAD 

• SCREEN SCALE 

• MISSILE RADAR BEAMWIDTH 

• RANGE GATE 

• MAXIMUM SHIP RADAR CROSS SECTION 

• MINIMUM SHIP RADAR CROSS SECTION 

• RADAR CROSS SECTION OF CHAFF 

• TRUE WIND DIRECTION 

• TRUE WIND SPEED 

• ALL PARAMETERS 

• RETURN TO MAIN PROGRAM 



4. Option Four 

This option terminates the program. 

E. SUBROUTINES 

A number of subroutines are made use of to break the program 
into smaller modules. Their respective functions and broad outlines 
are presented here to give an overview to the reader. 
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1 



Auto Sxibroutines 



There are two auto subroutines. The function of the first 
auto subroutine is to calculate the chaff launch bearing, BRGC, the 
first leg of the ship maneuvers, CSE2 and the direction to eneO^le 
the shorter turn from the initial course. In order to optimize miss 
distance, a fairly accurate estimate of the chaff launch bearing 
is obtained using the pre-computation method. The decision related 
to CSE2 calculation is based on the algorithm in section . The 
second subroutine works out the second leg of the maneuver and the 
direction to turn away from the chaff. 

2. Ship, Rescell And Chaff Subroutines 

These subroutines calculate the coordinates of the 

4 

corners of the resolution cell, the hexagon representing the chaff 
cloud and the polygon representing the ship. 

3. Shipcell And Chaffcell Subroutines 

These subroutines calculate the radar cross sections of 
the ship and the chaff inside the resolution cell to enable 
calculation of the centroid in the main program. 

4 . Aspect Subroutine 

The aspect angle of the ship as it exits the resolution 
cell is passed on to the Shipcell subroutine for subsequent RCS 
calculations. 

5. Missile, Wind, Path And Tadsle 1 Subroutines 

These are PLOT-88 subroutines. The first three are used 
to draw the missile line, wind line and paths traced by the 
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missile, ship and chaff. Table 1 sxibroutine constructs a table of 
essential parameters when Screnl & 2 subroutines are used. 

6. Scren 1 2 & 3 Subroutines 

Scren 1 subroutine draws the missile line of flight, wind 
vector and the ship polygon based on the initial conditions input 
by the operator. Scren 2 subroutine is used to draw the same 
picture but at 5 second intervals to offer a look at the situation 
dynamically. Scren 3 subroutine enables the operator to change 
desired parameters and updates the memory with their latest 
version. 
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V. RESULTS/RECOMMENDATIONS 



A. RESULTS 

The random variability of the parameters associated with the 
ship, missile and the chaff in a given battle scenario can result 
in a very large number of possible situations at sea. Some 
combinations of the parameters mentioned can produce a more 
favorable situation for the ship than others to thwart the anti- 
ship missile. In some cases, a single chaff round deployed by 
itself would have reduced effectiveness or might not be effective 
at all. A strong wind blowing from the direction of the missile or 
from the opposite side will move the chaff via the shortest path 
out of the resolution cell. In such cases, the chaff will exit the 
cell before the ship can do so; thus foiling the effort to break 
lock. Firing the chaff into the predicted range gate would avoid 
this problem. It is not considered appropriate to discuss results 
produced by the algorithm used in this work for all the possible 
cases. Some random situations have, therefore, been selected to 
give the reader an idea of the working of the model developed. 

1. Situation 1 

The salient input variables for this situation are as 

follows : 

• Missile Bearing - 000 Degrees 

• Missile Range - 10,000 Yards 

• Ship's Course - 225 Degrees 

• Wind Direction - 330 Degrees 

• Wind Speed - 15 Knots 
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The initial picture is shown in Figure 43. The bearing to launch 
as calculated by the program when selected in automatic mode is 
090. This has been worked out by the pre-computation technique. If 
the manual mode is selected instead, the user is presented the same 
picture without the launch bearing and asked to input the 
following: 

• Bearing to launch the chaff 

• First course to steer to exit the range gate 

• Direction to turn to first course 

• Second course to steer to get away from the missile line of 
fire 

• Direction to turn to second course 

In the automatic mode, however, the program works out the 
solution itself and then the result of the simulation can be viewed 
in two ways. The first graphic display shows the paths traced by 
the ship, the missile and the chaff. This is shown in Figure 44. 
In the case shown, the missile has missed the ship by 882.9 yards. 
The ship first maneuvered to steer southwards to break lock and 
then westwards to maximize miss distance. The second mode of 
presentation displays the situation every 5 seconds until the 
missile hits one of the two possible targets. The situation after 
15 seconds is shown in Figure 45. 
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Figure 43. Initial Conditions (Situation 1) 
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Figure 44. Paths display (Situation l) 
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Figure 45. Instantaneous Positions (Situation 1) 
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2 



Situation 2 



The input variables are as follows: 

• Missile Bearing - 280 Degrees 

• Missile Range - 10,000 Yards 

• Ship's Course - 150 Degrees 

• Wind Direction - 220 Degrees 

• Wind Speed - 15 Knots 

Figure 46 shows the initial conditions. The bearing to launch the 
chaff is 010. The paths traversed and the situation after 15 
seconds are shown in Figures 47 and 48 respectively. 
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Figure 46. Initial Conditions (Situation 2) 
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Figure 47. Paths display (Situation 2) 
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MINIMUM RANGE TO SHIP 906.4 TRRDS 
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Figure 48. Instantaneous Positions (Situation 2) 
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B. RECOMMENDATIONS 

The deployment of chaff in seduction mode poses a complicated 
problem in terms of the randomness of those variables that are not 
controllable by the ship. They include the missile's bearing, the 
beam width of its radar, the ship's course at the time of initial 
detection and the wind vector. Combined with the initial detection 
range of the missile these factors determine the size of the radar 
resolution cell and as a result, the time spent by the chaff round 
inside the cell. This is also the time duration available to the 
ship to vacate the cell and is, therefore, critical in solving the 
problem. 

The variables that can be controlled by the ship are the chaff 
range and bearing, ship's course and speed, its turn rates and 
acceleration. The last three have a limited contribution in this 
context. The control of the ship's heading and the bearing to fire 
the chaff is of the utmost significance. In order to increase the 
effectiveness of the control variables and to make the algorithm 
developed more adaptable, some recommendations are made in the 
ensuing paragraphs. 

1. Multiple Chaff Firing 

In order to tackle the cases where the chaff exits the 
radar resolution cell before the ship, multiple chaff rounds must 
be launched. The second chaff should be reseeded inside the cell 
prior to the exit of the first. Its initial placement should depend 
on how far the ship is from the cell bounds when the first round 
is about to exit or how much time it requires at the point to leave 
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the cell keeping the bloom time delay in mind. Thereafter, 
subsequent chaff rounds may be fired in the same way as required 
until the ship breaks lock. 

2. Chaff Walk Off 

This is another philosophy involving multiple chaff 
launches. If multiple chaff rounds are fired on the seune bearing 
but at increasing ranges, rapid break lock could be achieved and 
then the ship would need to carry out only one maneuver to increase , 
miss distance. The chaff rounds can be fired simultaneously or in 
quick succession so that they are waiting in anticipation and the 

I 

centroid shifts in a discrete manner. To explain this, let us study I 
an example shown in Figure 49. 

With the missile approaching from the north, the first 
chaff is placed at the cell's angle gate boundary on a bearing 
orthogonal to missile direction. The bearing can be offset to cater 
for the bloom time delay. If the chaff radar cross section is equal i 
to that of the ship, the centroid will be established at 
approximately the mid-point of the two, called XI. The resolution 
cell now moves to the right and is called Cell 'B'. At this point 
it includes the ship. Chaff 1 and Chaff 2. The centroid now shifts 
to X2 and the cell to Cell 'C'. The third chaff round displaces 
the cell position to Cell 'D' and the centroid to X3. It can be 
seen that the ship is now very close to the left hand angle gate 
boundary and can steer out of the cell to achieve a rapid break and 

I 

then increase the miss distance. We feel that this technique can 
enhance the effectiveness of chaff in anti-ship missile defense. 



78 



I 

MISSILE 

◄ ► CELL A 

CELLS 

CELLC 



CELLO 




Figure 49. Chaff Walk Off 

3. Multiple Solutions 

In order to achieve optimum results, the program can be 
used in manual mode to inject different control inputs. The program 
run time is quite short and a slight modification to compare 
results of both automatic and manual modes can recommend the best 
solution to be adopted in a particular situation. 

4. Battle Group Operation 

When operating in a battle group, non-judicious 
deployment of chaff can result in seducing the chaff onto a 
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friendly unit. It is, therefore, imperative to introduce safe 
firing bearings in the algorithm to ensure safety of ships in 
company . 

5. Fixed Chaff Launcher 

The program needs to be adapted to include units equipped 
with fixed chaff launchers to optimize firing bearings while 
minimizing maneuvers by the ship. The complexity of the problem is, 
of course, increased and the ship's actions are further limited. 
If there are limitations of CIWS firing bearings, they can also be 
included to allow for an integrated response. 

6. Integration with Active ECM 

If the ship is equipped with an active ECM capability, 
integration of chaff deployment with the active ECM response to an 
anti-ship missile threat can be a decisive factor. The break lock 
may first be achieved with active ECM and chaff can then be placed 
between the ship and the false target to deceive the missile and 
also to increase miss distance. Active decoys such as Nulka and 
Winnin employ similar techniques. 

7. Stochastic Modelling 

Due to the randomness of the variables involved in the 
anti-ship defense problem, stochastic modelling may be employed to 
achieve optimum results. In this paper, wind direction and speed 
are assumed to be constant. This is usually not true in real sea 
conditions. Noise in the wind vector can also be introduced through 
an appropriate probability distribution to depict more realistic 
environmental conditions. 
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8. Radar Cross Section Calculations 

The mathematical modelling to calculate radar cross 
sections of both the ship and the chaff, use some approximations 
when a simultaneous exit from the angle and range gates takes 
place. Also, in case of simultaneous spill-over from opposing sides 
of the cell through the range gate, the partial radar cross section 
is calculated assuming that the spill-over is from one side only. 
The cross section which lies outside the cell on the opposite side 
is thus neglected. The reader should be aware of these 
approximations and make suitable changes if so desired. 

9. Progreimming Modifications 

The simulation program is written in FORTRAN which has 
some inherent limitations. It is recommended that a more powerful 
and flexible language be used to improve upon the program 
efficiency and to incorporate the changes suggested more 
effectively. The dimension figure of the arrays can be increased 
to include larger initial detection ranges which may be applicable. 
Parameters such as chaff radius, ship's length, turning rates and 
acceleration/deceleration tables may be included as input variables 
to adapt the program to given chaff and ship charateristics . 

C. CONCLUSIONS 

Resurgence of decoys in recent years has highlighted their 
importance in the area of anti-ship missile defense. Although the 
missile radars are becoming smarter by including advanced ECCM 
features, chaff is considered to be a viable defensive technique 
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even if deployed in a stand alone manner such as centroid mode. 
This is especially valid for the smaller navies who have limited 
resources. Improvement in chaff deployment algorithms and, when 
possible, integrated operations with active ECM and CIWS assets can 
make chaff a potent defensive measure. It is simple to use and 
relatively inexpensive among the known type of decoys. Its 
intelligent use can make it effective and dependable. Although wars 
have not been won by a single factor alone, chaff can certainly 
become one of the important factors in a limited conventional 
conflict . 



82 



ooooooooooooooooooooo 



APPENDIX. SIMULATION PROGRAM 
******************************************* 



* DEPLOYMENT OF CHAFF IN CENTROID * 

* MODE AGAINST ANTI SHIP MISSILES * 

* USING A * 

* VARIABLE AZIMUTH AND ELEVATION * 

* LAUNCHER * 

* 24 AUG 89 SIDDIQ.FOR * 

* M. SIDDIQ * 



******************************************* 



THE PROGRAM HAS NO PRINTER/ PLOTTER OPTIONS 
ITS OUTPUT GOES TO THE MONITOR ONLY 

IT HAS AN AUTOMATIC OR MANUAL CHAFF LAUNCHING OPTION 
IT CREATES NO DATA FILE 

VARIABLES DECLARATION 



COMMON/TIME/ ITTG , TTG 
COMMON/TME/T 

COMMON/TRACK/ XM ( 60) , YM ( 60) , XS (60) , YS ( 60 ) , XC ( 60 ) , YC ( 60 ) 
COMMON/ WD/ GW, VW 
COMMON/MISS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) , R ( 60) , SC , VS ( 60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/EQ/ XET(60) , YET(60) 

COMMON/ SHCL/ AGS , RGS , BRGMS (60), S AA , PARTS , SAAN (60) 

COMMON/ CFCL/PARTC , CRAD, RGC , AGC 
COMMON/SHP/SL, SB 

COMMON/RC/RCSS , RCSSP ( 60 ) , RCSC , RCSCP ( 60 ) 

COMMON/MINN/RMIN 
COMMON/ CHF/BRGC, ISHOOT 
COMMON/TURN/TR, CSEOFF 
COMMON/FG/IFLAG1(60) ,IFLAG2(60) 

COMMON/ PRN/ I OPORT , MODEL 

COMMON/DOTS/XMD(60) ,YMD(60) ,XMDD(60) ,YMDD(60) ,XSD(60) ,YSD(60) , 
*XSDD(60) , YSDD(60) 

COMMON/REL/XR(60) ,YR(60) ,XRD(60) ,YRD(60) ,XRP(60) ,YRP(60) ,XRPD(60) , 
*YRPD(60) 

REAL AMM,NAV,T 
REAL CDIST 
REAL TTG,GT,X,Y 
REAL SIG,SIGD,AM,GM 
lOPORT =99 
MODEL =99 
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LIST OF VARIABLES 



C 

C 

C 

C 

C 

C 

C 



MISSILE: 

BRGM 

BRGMS 

BRGMT 



MISSILE (REF) 



C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



BRGMC 

R/R(N) 

RANGE 

R(1)/RMIN 

RMT(N) 

RMC 

VMM 

VM 

GM 

AM 

AMM 

XM 

YM 

XMD 

YMD 

XMDD 

YMD 

SHIP: 



SL 

SB 

CSE 

CSE2 

CSE3 

CSEOPT 

CSEDIFl 

TOPT 

GS 

VS 

AS (VR) 

TR 

U(ALSO U1,U2 

SAA 

XS 

YS 

XSD 

YSD 

XSDD 

YSSD 

RCSS 

RCSSMAX 

RCSSMIN 

RCSSP(l) 

RCSSP(N) 



BEARING OF MISSILE FROM SHIP(WRT NORTH) 
BEARING OF SHIP FROM MISSILE (WRT REFERENCE) 
BEARING OF EQUIVALENT TARGET FROM 

BEARING OF CHAFF FROM MISSILE (REF) 

INSTANTANEOUS MISSILE RANGE FROM SHIP 

INITIAL MISSILE RANGE 

SCALED INITIAL MISSILE RANGE 

MISSILE RANGE FROM EQUIVALENT TARGET 

MISSILE RANGE FROM CHAFF 

MISSILE VELOCITY 

SCALED MISSILE VELOCITY 

MISSILE DIRECTION 

MISSILE ACCELERATION 

SCALED MISSILE ACCELERATION (30/SC) 

X-COORD OF MISSILE POSITION 
Y-COORD OF MISSILE POSITION 
X-COORD OF MISSILE VELOCITY 
Y-COORD OF MISSILE VELOCITY 
X-COORD OF MISSILE ACCELERATION 
Y-COORD OF MISSILE ACCELERATION 



SHIP LENGTH 
SHIP BREADTH 
COURSE 

FIRST LEG OF NEW COURSE TO STEER 
SECOND LEG OF NEW COURSE TO STEER 
COURSE OPTION FOR FIRST/ SECOND LEG 
DIFFERENCE BETWEEN CSEl AND CSE2 
TIME TO TURN TO CSEOPT 
ANGLE BETWEEN REF AND SHIP'S HEAD 
SCALED SHIP VELOCITY 
SHIP ACCELERATION 
TURN RATE OF SHIP 
,U3)TURN DIRECTION OF SHIP 

ASPECT ANGLE OF SHIP (BETWEEN CSE & BETA) 

X-COORD OF SHIP POSITION 

Y-COORD OF SHIP POSITION 

X-COORD OF SHIP VELOCITY 

Y-COORD OF SHIP VELOCITY 

X-COORD OF SHIP ACCELERATION 

Y-COORD OF SHIP ACCELERATION 

RADAR CROSS SECTION OF SHIP 

MAXIMUM RADAR CROSS SECTION OF SHIP 

MINIMUM RADAR CROSS SECTION OF SHIP 

INITIAL RADAR CROSS SECTION OF SHIP 

INSTANTANEOUS RADAR CROSS SECTION OF SHIP 
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ANGLE BETWEEN REF USED FOR MICROSOFT 



C BETA 

FORTRAN & 

C 

C RESOLUTION CELL; 

C 

C RGATE 

C RG 

C AG 

C RGS 

CENTER 
C 
C 



C AGS 

CENTER 
C 
C 



BRGM 



RANGE GATE 
SCALED RANGE GATE 
ANGLE GATE 

DISTANCE BETWEEN RESCELL CENTER & SHIP 

ALONG THE RANGE GATE 
RGC 

CHAFF 

DISTANCE BETWEEN RESCELL CENTER & SHIP 

ALONG ANGLE GATE 
AGC 



C S2r/S2a 

CENTER 

C 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

CHAFF 

C 

C 

C 

C 

C 

C 

C 

C 

C 

OR 

C 

C 

C 

C 

TARGET 

C 

C 



Slr/Sla 

PARTS 

PARTC 

CHAFF: 



BRGC 

CDIST 

CRAD 

VC 

GC 

BRGC1&2 

RCSC 

RCSSC(l) 

RCSSC(N) 

XC 

YC 

EQUIVALENT TARGET 



BRGTM 

GET 



XET 

YET 

VET 



WIND: 



CHAFF 

DISTANCE BETWEEN RG/AG BOUNDARY & SHIP 

Rlr/Rla 

CHAFF 

PORTION OF HALF OF THE LONGITUDINAL AXIS 
INSIDE THE CELL 

PORTION OF THE SHIP INSIDE THE RESCELL 
CHAFF 



BEARING TO LAUNCH CHAFF 
CHAFF LAUNCH DISTANCE 
CHAFF RADIUS 
CHAFF SPEED 

DIRECTION OF MOTION OF CHAFF 

BEARINGS PERPENDICULAR TO BRGM TO FIRE 

RADAR CROSS SECTION OF CHAFF 
INITIAL RADAR CROSS SECTION OF CHAFF (=0) 
INSTANTANEOUS RADAR CROSS SECTION OF CHAFF 
X-COORD OF CHAFF POSITION 
Y-COORD OF CHAFF POSITION 



BEARING OF MISSILE FROM EQUIVALENT TARGET 
DIRECTION OF MOTION OF EQUIVALENT TARGET 

EQUIVALENT TARGET ANGLE 
X-COORD OF EQUIVALENT TARGET POSITION 
Y-COORD OF EQUIVALENT TARGET POSITION 
RATE OF CHANGE OR VELOCITY OF EQUIVALENT 

POSITION 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



GW WIND DIRECTION (FROM WHERE IT COMES FROM) 

VW TRUE WIND SPEED 

VWl SCALED WIND SPEED 

TW TRUE WIND DIRECTION (TO WHERE IT GOES) 

RELATIVE COORDINATES 



XR 

YR 

XRP 



YRP 



XRPD(l) 

YRPD(l) 

YRPD(N) 



VELOCITY 

C 

VELOCITY 

C 

DIRECTION 

C 

DIRECTION 

C 

DIRECTION 

C 

DIRECTION 



XRPD(N) 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 



XRPDl 

YRPDl 

XRPD2 

YRPD2 

XRPDD 

YRPDD 

MISCELLANEOUS; 

SC 

TTG 

ITTG 

XCENT 

YCENT 

NAV 

SIG 



SIGHT (CLOCKWISE) 



C 

C 

10 



SIGD 

PI = 4.0 * 
NAV =3.0 
T=0 

H = 1.0 
XCENT =4.0 



SHIP 

SHIP 

SHIP 



X-COORD OF RELATIVE POSITION OF 
Y-COORD OF RELATIVE POSITION OF 
X-COORD OF RELATIVE POSITION OF 
WITH THE RELATIVE COORDINATE 
SYSTEM ROTATED TO MISSILE BEARING 
Y-COORD OF RELATIVE POSITION OF SHIP 
WITH THE RELATIVE COORDINATE 
SYSTEM ROTATED TO MISSILE BEARING 
INITIAL X-COORD OF RELATIVE SHIP'S VELOCITY 
INITIAL Y-COORD OF RELATIVE SHIP'S VELOCITY 
INSTANTANEOUS Y-COORD OF REL SHIP'S 

INSTANTANEOUS X-COORD OF REL SHIP'S 



XRPl 


DISTANCE 


TRAVELLED 


BY 


SHIP 


IN 


POS 


XRP 


YRPl 


DISTANCE 


TRAVELLED 


BY 


SHIP 


IN 


POS 


YRP 


XRP2 


DISTANCE 


TRAVELLED 


BY 


SHIP 


IN 


NEG 


XRP 


YRP2 


DISTANCE 


TRAVELLED 


BY 


SHIP 


IN 


NEG 


YRP 



RELATIVE VELOCITY OF 
RELATIVE VELOCITY OF 
RELATIVE VELOCITY OF 
RELATIVE VELOCITY OF 
X-COORD OF RELATIVE 
Y-COORD OF RELATIVE 



SHIP IN XRPl DIRECTION 
SHIP IN YRPl DIRECTION 
SHIP IN XRP2 DIRECTION 
SHIP IN XRP2 DIRECTION 
ACCELERATION OF SHIP 
ACCELERATION OF SHIP 



SCREEN SCALE 
TIME TO GO 
PROGRAM RUN TIME 
X-COORD OF SCREEN CENTER 
Y-COORD OF SCREEN CENTER 
PROPORTIONAL NAVIGATION CONSTANT 
ANGLE BETWEEN REF AND 

RATE OF CHANGE OF LOS 



LINE OF 



ATAN(l.O) 
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YCENT =3.0 
SL = 200.0 
SB = 20.0 
IFLAGl(l) = 0 
IFLAG2(1) = 0 
IFLAG3 = 0 
C 

WRITE (*,*) 'INPUT MISSILE BEARING (000-359 DEGREES TRUE) ' 
READ(*,*) BRGM(l) 

WRITE(*,*) 'INPUT INITIAL MISSILE RANGE (< 15000 YARDS) ' 
READ(*,*) RANGE 

WRITE(*,*) 'INPUT SHIP HEADING (000-359 DEGREES TRUE)' 
READ(*,*) CSE(l) 

C MISSILE PARAMETERS 
C 

WRITE (*,*) 'INPUT MISSILE RADAR BEAMWIDTH (+/“ 1.0 TO +/“5.0 
♦DEGREES) ' 

READ(*,*) THETA 

WRITE(*,*) 'INPUT THE RANGE GATE (+/-10 TO +/-100 YARDS)' 
READ ( * , * ) RGATE 
C RADAR CROSS SECTION 
C 

C SHIP RCS ASSUMED MAX RIGHT AHEAD/ ASTERN & ABEAM AND MIN AT 
BOWS/QTRS 

WRITE (*,*)' INPUT RADAR CROSS SECTION OF SHIP & CHAFF' 

WRITE (*,*) ' FOR CHAFF TO BE EFFECTIVE THERE MUST BE AT LEAST' 
WRITE(*,*)' A 1:2 RCS RATIO ( CHAFF: SHIP) ' 

WRITE(*,*)' ' 

WRITE (*,*) 'INPUT MAX RADAR CROSS SECTION OF SHIP ' 

WRITE(*, *) '10000 - 100000 SQUARE METERS' 

WRITE(*,*)' ' 

READ(*,*) RCSSMAX 

WRITE (*,*) 'INPUT MINIMUM RADAR CROSS SECTION OF SHIP' 

WRITE (*,*) '5000 - 50000 SQUARE METERS' 

WRITE (*,*)' ' 

READ(*,*) RCSSMIN 

WRITE (*,*) 'INPUT RADAR CROSS SECTION OF CHAFF' 

WRITE (*,*) '5000 - 20000 SQUARE METERS ' 

WRITE (*,*)' ' 

READ(*,*) RCSC 
C WIND PARAMETERS 
C 

WRITE (*,*) 'INPUT WIND DIRECTION (000-359 DEGREES TRUE) ' 
READ(*,*) GW 

WRITE (*,*) 'INPUT TRUE WIND SPEED (0.0-50.0 KNOTS) ' 
READ(*,*) VW 

WRITE (*,*) 'INPUT SCREEN SCALE (YDS:1 INCH) ' 

WRITE(*,*)' RECOMMEND 200 YARDS TO THE INCH' 

READ(*,*) SC 
20 CONTINUE 
C SCALING 
C 
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C vs = SHIP'S SCALED VELOCITY CONVERTED TO YDS/SEC 

VS(1) = 20./(SC*1.8) 

C RG = SCALED RANGE GATE 

RG = RGATE/SC 

C VM = SCALED MISSILE VELOCITY 

VMM * 600. 

VM = VMM/(SC*1.8) 

C VWl = SCALED WIND VELOCITY 

VWl = VW/(SC*1.8) 

C SCALE THE RANGE FOR PLOTTING 

R(l) = RANGE/SC 
RMIN = R(l) 

C TR IS THE TURN RATE IN DEGREES/SEC 

TR = 5.0 

C VR IS THE ACCELERATION IN KNOTS/SEC*SC 

VR = 0.5/(SC*1.8) 

RCSS=RCSSMIN+ABS ( (RCSSMAX-RCSSMIN) *COS (2* (CSE(l) -BRGM(l) ) *PI/180) ) 
RCSSP(l) = RCSS 
RCSCP(l) = 0.0 
C TTG = TIME TO GO 

C 

TTG = RANGE* 1.8/VMM 

ITTG = IFIX(TTG +10*500 ./ (VM*SC) ) 

C MISSILE ACCELERATION 

C 

AMM=30./SC 

C SIG REF. TO HORIZONTAL (EAST) 

C 

SIG = (270. - BRGM(l) ) *PI/180. 

IF (BRGM(l) .GT.180. ) THEN 
GM = BRGM(l) - 180. 

ELSE 

GM = BRGM(l) + 180. 

END IF 

C ******INITIAL CONDITIONS ****** 

C MISSILE I. C. 

C 

XM(1) = XCENT + R(l) *SIN(BRGM(1) *PI/180. ) 

YM(1) = YCENT + R(l) *COS(BRGM(l) *PI/180. ) 

XMD(l) = VM*SIN(GM* PI/180.) 

YMD(l) = VM*COS(GM* PI/180.) 

C SHIP I. C. 

C 

XS(1) = XCENT 
YS(1) = YCENT 

XSD(l) = VS(1) * SIN(CSE(1) *PI/180. ) 

YSD(l) = VS(1) * COS(CSE(1)*PI/180. ) 

XSDD(l) = 0.0 
YSDD(l) = 0.0 
GS = ATAN2 (YSD(l) ,XSD(1) ) 

C FIND THE NEW COURSE TO STEER AND BEARING TO LAUNCH CHAFF 



88 



c 

30 WRITE (*,*)' INPUT (1) FOR AUTOMATIC CHAFF LAUNCH ' 
WRITE(*,*)' AND NEW SHIP COURSE ' 

WRITE(*,*)' ' 

WRITE (*,*)' INPUT (2) FOR MANUAL CHAFF LAUNCH ' 

WRITE (*,*)' AND NEW SHIP COURSE ' 

READ(*,*) ISHOOT 
IF (ISHOOT. EQ.l) THEN 
CALL AUTOl (CSE2,U2) 

U = U2 

CALL AUT02 (CSE3 ,U3 , CSE2 ) 

U1 = U3 

C DRAW THE INITAL MISSILE RESOLUTION CELL 

C 

WRITE(*,*) '**** MISSILE DETECTION **** ' 

RMT(l) = R(l) 

BRGTM(l) = BRGM(l) 

XET(l) = XS(1) 

YET(l) = YS(1) 

CALL SCRENl 

ELSE IF ( ISHOOT. EQ. 2) THEN 

WRITE(*, *) '**** MISSILE DETECTION **** ' 

RMT(l) = R(l) 

BRGTM(l) = BRGM(l) 

XET(l) = XS(1) 

YET(l) = YS(1) 

CALL SCRENl 

WRITE (*,*) 'INPUT TRUE BEARING TO LAUNCH CHAFF' 

READ(*,*) BRGC 

WRITE (*,*) 'INPUT FIRST LEG OF NEW SHIP COURSE TO STEER ' 
READ(*,*) CSE2 

40 WRITE(*,*)' INPUT (1) TO TURN THE SHIP RIGHT ' 

WRITE(*,*)' INPUT (-1) TO TURN THE SHIP LEFT ' 

READ(*,*) U2 
U = U2 

IF ( (U2.NE.1.0) .AND. (U2.NE.-1.0) ) THEN 
GO TO 40 
END IF 

WRITE (*,*) 'INPUT SECOND LEG OF NEW SHIP COURSE TO STEER ' 
READ(*,*) CSE3 

50 WRITE(*,*)' INPUT (1) TO TURN THE SHIP RIGHT ' 

WRITE(*,*)' INPUT (-1) TO TURN THE SHIP LEFT ' 

READ(*,*) U3 
U1 = U3 

IF ((U3.NE.1.0) .AND. (U3.NE.-1.0)) THEN 
GO TO 50 
END IF 
ELSE 

GO TO 30 
END IF 

C CHAFF INITIAL CONDITIONS 
C 
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on non ooo 



C CDIST = DISTANCE FM SHIP CHAFF IS LAUNCHED 
C CRAD = CHAFF RADIUS 

C CHAFF IS FIRED TO THE EDGE OF THE CELL BUT THE DISTANCE IS 
REDUCED 

C BY ONE CHAFF RADIUS TO CATER FOR THE SHIP'S MOVEMENT DURING THE 
C BLOOM TIME DELAY 

C VC = CHAFF SPEED IN YDS/SEC 

CRAD = 100. /SC 

CDIST = (R(l) *TAN( THETA * PI/180)) - CRAD 
VC = VW / (SC*1.8) 

GC = GW + 180. 

IF (GC.GT.360) THEN 
GC = GC - 360 
ENDIF 

XC(1) = XS(1) 

YC(1) = YS(1) 



************ main loop ****************** 

DO 100 N = 1,ITTG 

IF (IFLAGl (N-1) .EQ. 1) GOTO 110 



IFLAGl(N) = 0 
IFLAG2(N) = 0 

RANGES AND ANGLES (WRT HORIZONTAL REFERENCE 



R(N) = SQRT( (XS(N)-XM(N) ) **2+{YS(N)-YM(N) )**2) 

RCM = SQRT( (XC(N) -XM(N) ) **2+(YC(N) -YM(N) ) **2) 

RMT(N) = ((XM(N)-XET(N))**2+(YM(N)-YET(N))**2)**0.5 
IF( ( (XET(N) -XM(N) ) .EQ.O. ) .AND. { (YET(N) -YM(N) ) .EQ.O. ) ) THEN 
BRGMT(N) = 0.0 
ELSE 

BRGMT(N) = ATAN2 (YET(N) -YM(N) ,XET(N) -XM(N) ) 

END IF 

IF( ( (XET(N) -XM(N) ) .EQ.O.) .AND. ( ( YET (N) -YM (N) ) .EQ.0.))THEN 
BRGTM(N) = 0.0 
ELSE 

BRGTM(N) = (180./PI) *ATAN2 (XM(N) -XET(N) , YM(N) -YET(N) ) 
END IF 

IF(((XS(N)-XM(N)) .EQ.O.) .AND. ( ( YS (N) -YM (N) ) .EQ.O.)) THEN 
BRGMS(N) =0.0 
ELSE 

BRGMS(N) = ATAN2(YS(N)-YM(N) ,XS(N)-XM(N) ) 

END IF 

IF( { (XC(N) -XM(N) ) .EQ.O. ) .AND. ( (YC(N) -YM(N) ) .EQ.O. ) ) THEN 
BRGMC =0.0 
ELSE 

BRGMC = ATAN2 (YC(N) -YM(N) ,XC(N) -XM(N) ) 

END IF 

C RANGE AND ANGLE GATES 
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c 

AGS = R(N) * SIN(ABS(BRGMS(N)-BRGMT(N) ) ) 

AGC = RCM * SIN(ABS(BRGMC-BRGMT(N) ) ) 

RGS = ABS(R(N) * COS (ABS (BRGMS (N) -BRGMT (N) ) ) -RMT (N) ) 

RGC = ABS (RCM * COS (ABS (BRGMC-BRGMT(N) ) ) -RMT (N) ) 

AG = RMT(N) * TAN (THETA* PI/ 180.) 

C RCS MODEL 
C 

RCSS=RCSSMIN+ABS ( (RCSSMAX-RCSSMIN) *COS (2* (CSE(N) -BRGM(N) ) 

* *PI/180)) 

RCSSP(N) = RCSS 

IF ( (RGS.GT. (RG-SL/(SC*2.0) ) ) .OR. 

* (AGS.GT. (AG-SL/(SC*2.0) ) ) ) THEN 

CALL SHPCLL(N) 

RCSSP(N) = RCSS * PARTS 
END IF 

C CHAFF RCS IS ASSUMED TO BE ZERO UNTIL FULL BLOOM 

C 

IF (N.GT.5) THEN 
RCSCP(N) = RCSC 

IF ( (RGC.GT. (RG-CRAD) ) .OR. (AGC.GT. (AG-CRAD) ) ) THEN 
CALL CHFCLL 

RCSCP(N) = RCSC * PARTC 
END IF 
ELSE 

RCSCP(N) =0.0 
ENDIF 

C EQUIVALENT TARGET 
C 

IF ( (RCSSP(N) .EQ.0.0) .AND. (RCSCP(N) .EQ.0.0) ) THEN 
XET(N+1) = XET(N) 

YET(N+1) = YET(N) 

ELSE 

XET(N+1) = (RCSSP(N) *XS(N)+RCSCP(N) *XC(N) )/(RCSSP(N)+ 

* RCSCP(N)) 

YET(N+1) = (RCSSP(N) *YS (N)+RCSCP(N) *YC(N) )/ (RCSSP(N)+ 

* RCSCP(N)) 

END IF 

C ANGLE BETWEEN HORIZONTAL REFERENCE AND DIRECTION OF MOVEMENT 
OF 

C EQUIVALENT TARGET 

IF (N.GE.2) THEN 

IF( ( (XET(N) -XET(N-l) ) .EQ.O. ) .AND. ( (YET(N) -YET(N-l) ) .EQ. 0. ) ) 

* THEN 

GET =0.0 
ELSE 

GET = ATAN2 (YET(N)-YET(N-l) ,XET(N) -XET(N-l) ) 

END IF 

C EQUIVALENT TARGET VELOCITY 

C 

VET = SQRT((XET(N)-XET(N-1))**2+(YET(N)-YET(N-1))**2)/H 
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ELSE 

GET = GS 
VET = VS(N) 

END IF 

C MISSILE NAVIGATION DYNAMICS 

C 

IF( ( (XET(N)-XM(N) ) .EQ.O.) .AND. ( (YET(N) -YM(N) ) .EQ.O.) )THEN 
SIG =0.0 
ELSE 

SIG = ATAN2( (YET(N)-YM(N) ) , (XET(N)-XM(N) ) ) 

END IF 

IF (RMT(N) .EQ.0.0) THEN 
SIGD =0.0 
ELSE 

SIGD = (VET*SIN(GET-SIG) -VM*SIN(GM-SIG) )/RMT(N) 

END IF 

AM=SIGD*VM*NAV 
IF(AM.GT.AMM) THEN 
AM=AMM 

ELSE IF(AM.LT.-AMM) THEN 
AM = -AMM 
END IF 

C TO TURN OFF MISSILE SEEKER 

C 

IF (RMT(N) .LT. (250./SC) ) THEN 
AM = 0.0 

IF (RMT(N-l) .LT.RMT(N) ) THEN 
ITTG = N 
IFLAGl(N) = 1 
END IF 
END IF 

C MISSILE EQUATIONS 

C 

XMDD(N)= -AM*SIN(SIG) 

YMDD(N)= AM*COS(SIG) 

IF( (YMD(N) .EQ.O. ) .AND. (XMD(N) .EQ.O. )) THEN 
GM = 0.0 
ELSE 

GM = ATAN2 (YMD(N) ,XMD(N) ) 

END IF 

XM(N+1)= XM(N) + XMD(N)*H + 0 . 5*XMDD(N) *H**2 
YM(N+1)= YM(N) + YMD(N)*H + 0 . 5*YMDD (N) *H**2 
XMD(N+1)= VM*COS(GM) + XMDD(N)*H 
YMD(N+1)= VM*SIN(GM) + YMDD(N)*H 
C SHIP EQUATIONS 

C 

IF (BRGM(N) .LT.0.0) THEN 
BRGM(N) = BRGM(N) + 360. 

END IF 

XSDD(N) = -VS(N)*(PI*TR/180.)*U1*SIN(CSE(N)*PI/180.) + 

VR*SIN 

* (PI*CSE(N)/180. ) 
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YSDD(N) = VS(N)*(PI*TR/180.)*U1*COS(CSE(N)*PI/180.) + 

VR*COS 

* (PI*CSE(N)/180. ) 

XSD(N) = VS(N) *SIN(PI*CSE(N)/180.) 

YSD(N) = VS(N)*COS(PI*CSE(N)/180.) 

IF (VS (N) .LT. (25./(1-8*SC) ) ) THEN 
VS(N+1) = VS(N) + VR*H 
ELSE 

VS(N+1)=25./(1*8*SC) 

END IF 

XS(N+1) = XS(N) + XSD(N)*H + 0.5*XSDD(N) *H**2 
YS(N+1) = YS(N) + YSD(N)*H + 0. 5*YSDD(N) *H**2 
IF(((XM(N+1)-XS(N+1)) .EQ.O.) .AND. ( (YM(N+1) -YS (N+1) ) .EQ.O.)) 

* THEN 

BRGM(N+1) =0.0 
ELSE 

BRGM(N + 1) 

(180./PI)*ATAN2 (XM(N+1)-XS(N+1) , YM(N+1) -YS (N+1) ) 

END IF 

CHAFF DYANAMICS AND RELATIVE COORDINATES EQUATIONS 



IF (N.LE.5) THEN 

CHAFF'S ARBITRARY POSITION BEFORE IT BLOOMS 



XC(N+1) = XS(N+1) + CDIST*SIN(BRGC*PI/180.) 
YC(N+1) = YS(N+1) + CDIST*COS(BRGC*PI/180.) 

ELSE 

XC(N+1) = XC(N) + VC*SIN(GC*PI/180. ) 

YC(N+1) = YC(N) + VC*COS(GC*PI/180. ) 

ENDIF 

SHIFT OF THE ORIGIN OF COORDINATE SYSTEM TO THE CHAFF 



XR(N) = XS(N)-XC(N) 

YR(N) = YS(N)-YC(N) 

ROTATION OF THE COORDINATE SYSTEM TO THE MISSILE DIRECTION 



XRP(N) = (XR(N)*COS(BRGM(N)*PI/180)) - (YR(N)* 

* SIN(BRGM(N) *PI/180) ) 

YRP(N) = (XR(N) *SIN(BRGM(N) *PI/180) ) + (YR(N)* 

* COS(BRGM(N) *PI/180) ) 

T=T+H 

C ALTERATION OF COURSES 
C 

C WATCH FOR CSE(N) - CSE2/CSE3 NEVER BEING LESS THAN TR 

C WHILE THE SHIP IS WITHIN THE RG BOUNDS UPDATE COURSE TO ALTER 

TO 

C THE FIRST LEG UNTIL THE TURN IS COMPLETED 
IF (ABS(YRP(N) ) .LE.RG) THEN 

IF ( (ABS(CSE(N) -CSE2) ) .LT.TR) THEN 
U = 0.0 
END IF 

CSE(N+1) = CSE(N) + TR*H*U 
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C WHILE THE SHIP IS WITHIN THE RG BOUNDS BUT IT HAS EXITED THE 
AG 

C BOUNDS UPDATE COURSE TO ALTER TO THE FIRST LEG UNTIL THE TURN 
C IS COMPLETED 

IF (ABS(XRP(N) ) .GT.AG) THEN 

IF ( (ABS(CSE(N)-CSE3) ) .LT.TR) THEN 

U1 = 0.0 

ENDIF 

CSE(N+1) = CSE(N) + TR*H*U1 
ENDIF 
ELSE 

C IF THE SHIP IS OUT OF THE RG BOUNDS UPDATE COURSE TO ALTER TO 
THE 

C SECOND LEG UNTIL THE TURN IS COMPLETED 

IF ((ABS(CSE(N)-CSE3)) .LT.TR) THEN 
U1 = 0.0 
ENDIF 

CSE(N+1) = CSE(N) + TR*H*U1 
ENDIF 

IF (CSE(N+1) .GE.360. ) THEN 
CSE(N+1) = CSE(N+1) - 360.0 
END IF 

IF (CSE(N+1) .LT.O. ) THEN 
CSE(N+1) = 360.+CSE(N+l) 

END IF 

DETERMINE IF MISSILE HITS THE SHIP 



IF( (R(N)*SC) .LT.50.0)THEN 

WRITE (*,*) 'RANGE UNDER 50 YARDS MISSILE HITS SHIP' 
END IF 

RMIN = AMIN1(RMIN,R(N) ) 

IF (RMT(N) .LT. (500./SC) ) THEN 
IFLAG2(N) = 1 
H = 0.1 
END IF 



100 CONTINUE 
110 CONTINUE 

TO PREVENT ANY UNDEFINED VALUES IF ITTG <60 



DO 120 I=(ITTG+1) ,60 
XC(I) = XC(I-l) 
XC(I) = XC(I-l) 
YC(I) = YC(I-l) 
XS(I) = XS(I-l) 
YS(I) = YS(I-l) 
XM(I) = XM(I-l) 
YM(I) = YM(I-l) 
XET(I) = XET(I-l) 
YET(I) = YET(I-l) 
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CONTINUE 



C 
C 
C 

RMIN = SC*RMIN 

Qic ic it it ic it "k it it it ic ic ic ic it it it it it If it it it it it it it it irk ft it -k it ii it it ic it it it it it it ii ii ic irk it it it it "k it it It 

C PRINT & PLOT 

C 

130 WRITE (*,*) 'INPUT (1) 

WRITE(*,*)' ' 

WRITE(*,*) 'INPUT (2) 

WRITE(*,*) ' 

WRITE (*,*)' ' 

WRITE (*,*) 'INPUT (3) 

WRITE (*,*)' ' 

WRITE (*,*) 'INPUT (0) 

WRITE (*,*)' ' 

READ(*,*) IDRW 
IF(IDRW.EQ.l) THEN 
CALL PATH 
ELSE IF(IDRW.EQ.2) THEN 
CALL SCREN2 

ELSE IF (IDRW. EQ. 3) THEN 
CALL SCREENS (IFLAG3) 

IF (IFLAG3.EQ.1) THEN 
GO TO 10 
ENDIF 

IF (IFLAG3.EQ.2) THEN 
PI = 4.0 * ATAN(l.O) 

NAV =3.0 
T=0 

H = 1.0 
XCENT =4.0 
YCENT =3.0 
SL = 200.0 
SB = 20.0 
IFLAGl(l) = 0 
IFLAG2(1) = 0 
GO TO 20 
ENDIF 

ELSE IF(IDRW.EQ.O) THEN 
GO TO 999 
END IF 
GO TO 130 
999 STOP 
END 
C 

c* *************************************************************** 
SUBROUTINE AUTOl (CSE2,U2) 

C THIS SUBROUTINE CALCULATES THE BEARING TO LAUNCH THE CHAFF AND 
THE 

C FIRST LEG OF THE NEW COURSES TO MANEUVER THE SHIP 



TO DRAW MISSILE & SHIP & CHAFF PATH' 

TO DRAW RESOLUTION CELL IN ' 

FIVE SECOND INCREMENTS ' 

TO CHANGE SCENARIO PARAMETERS' 

TO QUIT ' 
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COMMON/TIME/ITTG , TTG 
COMMON/WD/GW , VW 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/MISS/RG , THETA , AG , PI 
COMMON/ CHF/BRGC, ISHOOT 

COMMON/REL/XR(60) ,YR(60) ,XRD(60) ,YRD(60) ,XRP(60) ,YRP(60) ,XRPD(60) , 
*YRPD(60) 

COMMON/TURN/TR, CSEOFF 

REAL CSEOPT(4) ,TOPT(4) ,CSE2 ,CSEDIF(4) ,CDIF1, ACC, VSMAX 
REALXRPl , YRPl , XRP2 , YRP2 , XRPDl , YRPDl , XRPD2 , YRPD2 , XRPDD, YRPDD 
INTEGER TTG1,T 
REAL CDIST 

POSSIBLE CHAFF LAUNCH BEARINGS 



BRGCl = BRGM(l) + 90 
BRGC2 = BRGM(l) - 90 
IF ( BRGCl. GT. 360) THEN 
BRGCl = BRGCl - 360 
ELSEIF (BRGCl. LT.O) THEN 
BRGCl = BRGCl +360 
ENDIF 

IF (BRGC2.GT.360) THEN 
BRGC2 = BRGC2 - 360 
ELSEIF (BRGC2.LT.0) THEN 
BRGC2 = BRGC2 +360 
ENDIF 

WRITE (*,*) 'BRGC1=' , BRGCl 
WRITE (*,*) 'BRGC2=' ,BRGC2 
GC = GW + 180 
VWl = VW/(SC*1.8) 

COURSE OPTIONS FOR THE FIRST LEG 



CSEOPT(l) = BRGM(l) 

CSEOPT(2) = BRGM(l) + 180 
COURSE OPTIONS FOR THE SECOND LEG 



CSEOPT(3) = BRGM(l) + 90 
CSEOPT(4) = BRGM(l) - 90 
DO 10 1=1,4 

IF (CSEOPT(I) .GT.360) THEN 
CSEOPT(I) = CSEOPT(I) - 360 
ELSEIF (CSEOPT(I) .LT.O) THEN 
CSEOPT(I) = CSEOPT(I) +360 
ENDIF 
CONTINUE 

CALCULATION OF INITIAL VELOCITY OF RELATIVE VELOCITY 



TEMP1=(VS(1)*(SIN( (CSE(l)-BRGM(l) )*PI/180) ) ) 
TEMP2=(VW1* (SIN( (GW-BRGM(l) ) +PI/180) ) ) 
XRPD(l) = TEMP1+TEMP2 
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TEMP3=(VS(1)*(COS((CSE(1)-BRGM(1))*PI/180) ) ) 

TEMP4=(VW1* (COS ( (GW-BRGM(l) ) *PI/180) ) ) 

YRPD(l) = TEMP3+TEMP4 

COURSE DIFFERENCE BETWEEN INITIAL COURSE AND THE COURSE OPTIONS 



DO 20 I = 1,4 

CSEDIF(I) = ABS(CSE(1) - CSEOPT(I)) 

IF (CSEDIF(I) .GT.180) THEN 
CSEDIF(I) = 360 - CSEDIF(I) 

ENDIF 

CONTINUE 

TIME TAKEN TO TURN TO COURSE OPTIONS 3 i 4 



TOPT(3) = CSEDIF(3)/TR 
TOPT(4) = CSEDIF(4)/TR 
INITIAL VALUES OF VARIABLES INVOLVED 



VSMAX = 25/(1.8*SC) 

TTGl = IFIX(TTG) 

ACC = VSMAX * TR * PI/ 180 
T = 1 

CRAD = 100. /SC 

CDIST = (R(l) * TAN(THETA*PI/180) ) - CRAD 
PRE-COMPUTATION TO CALCULATE CHAFF LAUNCH BEARING 



XRPl = CDIST 
YRPl =0.0 
XRP2 = -XRPl 
YRP2 =0.0 
XRPDl = XRPD(l) 

YRPDl = YRPD(l) 

XRPD2 = XRPD(l) 

YRPD2 = YRPD(l) 

XRPDD =0.0 
YRPDD =0.0 

IF COURSE IS ALTERED TOWARDS COURSE OPTION 3 (TO THE RIGHT) 



DO 30 I = 1,TTG1 

IF (I.LT.INT(TOPT(3) ) ) THEN 

XRPDD = ACC * COS(ABS(CSE(l) - BRGM(l)) * PI/180) 
YRPDD = -ACC * SIN(ABS(CSE(1) - BRGM(l)) * PI/180) 
ELSE 

XRPDD =0.0 
YRPDD = 0.0 
ENDIF 

XRPl = XRPl + XRPD1*T + 0.5 * XRPDD*T**2 
YRPl = YRPl + YRPD1*T + 0.5 * YRPDD*T**2 
XRPDl = XRPDl + XRPDD 
YRPDl = YRPDl + YRPDD 
30 CONTINUE 

IF COURSE IS ALTERED TOWARDS COURSE OPTION 4 (TO THE LEFT) 
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DO 40 I=1,TTG1 

IF (I.LT.INT(TOPT(4) ) ) THEN 

XRPDD = -ACC * COS(ABS(CSE(l) - BRGM(l) ) * PI/180) 

YRPDD = ACC * SIN(ABS(CSE(1) - BRGM(l)) * PI/180) 

ELSE 

XRPDD =0.0 
YRPDD =0.0 
ENDIF 

XRP2 = XRP2 + XRPD2*T + 0.5 * XRPDD*T**2 
YRP2 = YRP2 + YRPD2*T + 0.5 * YRPDD*T**2 
XRPD2 = XRPD2 + XRPDD 
YRPD2 = YRPD2 + YRPDD 
40 CONTINUE 

C CALCULTAION OF CHAFF LAUNCH BEARING 

C 

IF (ABS(XRPl) .GT.ABS(XRP2) ) THEN 
BRGC = BRGC2 
ELSE 

BRGC = BRGCl 
ENDIF 

C CALCULATION OF THE FIRST LEG OF THE NEW COURSES 

C 

IF (YRPD(l) .GT.O) THEN 
CSE2 = CSEOPT(l) 

ELSE 

CSE2 = CSEOPT(2) 

ENDIF 

WRITE (*,*) 'CSE2=' ,CSE2 

C COURSE DIFFERENCE BETWEEN INITIAL COURSE AND THE FIRST LEG 

C 

CDIFl = ABS(CSE(1) - CSE2) 

C CALCULATION OF TURN DIRECTION INDICATOR FOR THE FIRST LEG 

C 

IF (CDIFl. GE. 180) THEN 

IF (CSE2.GT.CSE(1) ) THEN 
U2 = -1.0 
ELSE 

U2 = 1.0 
ENDIF 

ELSEIF (CDIF1.lt. 180) THEN 
IF (CSE2.GT.CSE(1) ) THEN 
U2 = 1.0 
ELSE 

U2 = -1.0 
ENDIF 
ENDIF 
RETURN 
END 

C *************************************************************** 
SUBROUTINE AUT02 (CSE3 ,U3 , CSE2) 

C THIS SUBROUTINE CALCULATES THE SECOND LEG OF THE NEW COURSES 
C TO MANEUVER THE SHIP 
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COMMON/MISS/RG , THETA , AG , PI 
COMMON/ CHF/BRGC, ISHOOT 

COMMON/REL/XR(60) ,YR(60) ,XRD(60) ,YRD(60) ,XRP(60) ,YRP(60) ,XRDP(60) , 
*YRDP(60) 

COMMON/ABCD/BRGM(60) ,R(60) ,SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

REAL CSEOPT ( 4 ) , CSE3 , BRGM 
CSEOPT(l) = BRGM(l) 

CSEOPT(2) = BRGM(l) + 180 

CSEOPT(3) = BRGM(l) + 90 

CSEOPT(4) = BRGMU) " 90 

DO 10 1-1,4 

IF (CSEOPT(I) .GT.360) THEN 
CSEOPT(I) = CSEOPT(I) - 360 
ELSEIF (CSEOPT(I) .LT.O) THEN 
CSEOPT(I) = CSEOPT(I) + 360 
ENDIF 

10 CONTINUE 

BRGCl = BRGM(l) + 90 
BRGC2 = BRGM(l) - 90 
IF (BRGCl. GT.360) THEN 
BRGCl = BRGCl - 360 
ELSEIF (BRGCl. LT.O) THEN 
BRGCl = BRGCl +360 
ENDIF 

IF (BRGC2. GT.360) THEN 
BRGC2 = BRGC2 - 360 
ELSEIF (BRGC2.LT.0) THEN 
BRGC2 = BRGC2 + 360 
ENDIF 

CALCULATION OF THE SECOND LEG OF THE NEW COURSES 



IF (BRGC.EQ. BRGCl) THEN 
CSE3 = CSEOPT (4) 

ELSE 

CSE3 = CSEOPT (3) 

ENDIF 

WRITE (*,*)' CSE3= ' , CSE3 

C CALCULATION OF THE TURN DIRECTION INDICATOR FOR 
C THE SECOND LEG OF THE NEW COURSES 

C 

IF ( CSE3.EQ. CSEOPT ( 3) ) THEN 
IF (CSE2. EQ. CSEOPT ( 1) ) THEN 
U3 = 1.0 

ELSEIF ( CSE2.EQ. CSEOPT (2 ) ) THEN 
U3 = -1.0 
ENDIF 

ELSEIF (CSE3. EQ. CSEOPT (4 ) ) THEN 
IF ( CSE2.EQ. CSEOPT ( 1) ) THEN 
U3 = -1.0 

ELSEIF (CSE2. EQ. CSEOPT (2 ) ) THEN 
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1.0 



U3 = 

ENDIF 
ENDIF 
RETURN 
END 

******** 

SUBROUTINE PATH 

C THIS SUBROUTINE OF PLOT-88 DISPLAYS THE PATHS FOLLOWED 
C BY THE SHIP, MISSILE AND CHAFF 
COMMON/ PRN/ lOPORT , MODEL 
COMMON/TIME/ ITTG , TTG 

COMMON/TRACK/ XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 
COMMON/MINN/RMIN 
REAL XCP(55) ,YCP(55) 

CHARACTER*20 LS,U1,LC,YD 
CHARACTER* 2 5 LRMIN 
YD = 'YARDS' 

LRMIN=' MINIMUM RANGE TO SHIP' 

SHIFT OF CHAFF POSITION ARRAY TO PLOT CHAFF AFTER IT BLOOMS 



DO 10 I = l,(ITTG-5) 

XCP(I) = XC(I+5) 

YCP(I) = YC(I+5) 

CONTINUE 

TO PREVENT UNDEFINED VALUES IF ITTG-5 < 55 



DO 20 I=(ITTG-4) ,55 
XCP(I) = XCP(I-l) 

YCP(I) = YCP(I-l) 

20 CONTINUE 
LS = 'S' 

LM = 'M' 

LC = 'C' 

K = 2*ITTG/3 

CALL PLOTS (0,IOPORT, MODEL) 

CALL CURVE (XM,YM, -ITTG, 0. 1) 

CALL SYMBOL(XM(K) ,YM(K) , 0 . 20 , LM, 0 . 0 , 1) 

CALL CURVE (XS,YS, -ITTG, 0. 1) 

CALL SYMBOL(XS(l) ,YS(1) , 0 . 20 , LS , 0 . 0 , 1) 

CALL CURVE (XCP,YCP, -ITTG, 0.1) 

CALL SYMBOL(XC(6) ,YC(6),0.20,LC,0.0,1) 

CALL SYMBOL(0.5,0.2,0.25,LRMIN,0.0,21) 

CALL NUMBER(5.5,0.2, 0.25,RMIN,0.0, 1) 

CALL SYMBOL(7.0,0.2, 0.25,YD,0.0,5) 

CALL PLOT(0, 0,999) 

RETURN 

END 

Q* ************************************************************** * 

SUBROUTINE SCRENl 

C THIS PLOT-88 SUBROUTINE DISPLAYS THE INITIAL CONDITIONS 
C OF THE SHIP AND THE MISSILE 
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COMMON/ PRN/ lOPORT , MODEL 
COMMON/ CHF/BRGC, ISHOOT 
CHARACTER* 2 5 CL 

CL = 'AUTO CHAFF LAUNCH BEARING' 

CALL PLOTS (0,IOPORT, MODEL) 

CALL SHIP(l) 

CALL MISSLE(l) 

CALL RESCLL(l) 

CALL WIND 
CALL TABLEl(l) 

IF (ISHOOT.EQ.l) THEN 

CALL SYMBOL(1.0,5.7,0.16,CL,0.0,25) 

CALL NUMBER (6.0,5.7,0.16, BRGC , 0 . 0 , - 1 ) 

END IF 

CALL PLOT(0, 0,999) 

RETURN 

END 

Ct<^*************************************************************** 
SUBROUTINE SCREN2 

C THIS PLOT-88 SUBROUTINE DISPLAYS THE INSTANTANEOUS POSITIONS 
C OF THE SHIP, CHAFF AND THE MISSILE EVERY 5 SECONDS 
COMMON/PRN/IOPORT, MODEL 
COMMON/TIME/ ITTG, TTG 
COMMON/FG/IFLAG1(60) ,IFLAG2(60) 

L = ITTG/ 5 
DO 31 J=1,L 
K = J * 5 

IF (IFLAG2(N) .EQ.l) THEN 
K = J+1 
END IF 

CALL PLOTS ( 0 , lOPORT , MODEL) 

CALL SHIP(K) 

CALL MISSLE(K) 

IF (IFLAGl(K) .EQ.O) THEN 
CALL RESCLL(K) 

END IF 
CALL WIND 

C CHAFF NOT DISPLAYED FOR THE FIRST 5 SECONDS TO SHOW BLOOM TIME 

DELAY 

C 



IF (J.GT.l) THEN 
CALL CHAFF (K) 

ENDIF 

CALL TABLEl(K) 

CALL PLOT (0, 0,999) 

31 CONTINUE 
RETURN 
END 

C** **************************************************** ********** 
SUBROUTINE SHIP(J) 
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C THIS SUBROUTINE CALCULATES THE COORDINATES OF THE POLYGON 
C REPRESENTING THE SHIP 

COMMON/TRACK/ XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 
COMMON/MIS S/RG , THETA, AG , PI 

COMMON/ABCD/BRGM(60) ,R(60) , SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/SHP/SL, SB 
REAL XSHIP(5) ,YSHIP(5) 

C ROTATES THE SHIP DEPENDENT UPON THE COURSE 

C 

XSHIP(1)=XS(J) + (SL/2.)*SIN(CSE(J)*PI/180.)/SC 
XSHIP(2)=XS(J) + (SL/4.)*SIN(CSE(J)*PI/180.)/SC + 

* (SB/2 . ) *SIN ( (CSE(J) +90. ) +PI/180. )/SC 
XSHIP(3)=XS(J) - (SL/2.)*SIN(CSE(J)*PI/180.)/SC + 
*(SB/2.)*SIN( (CSE(J)+90.)*PI/180. )/SC 
XSHIP(4)=XS(J) - (SL/2.)*SIN(CSE(J)*PI/180.)/SC + 

* (SB/2 . ) *SIN ( (CSE ( J) -90 . ) *PI/180 . ) /SC 
XSHIP(5)=XS(J) + (SL/4. ) *SIN(CSE(J) *PI/180. )/SC + 

* (SB/2. )*SIN( (CSE(J)-90.) *PI/180. )/SC 
YSHIP(1)=YS(J) + (SL/2.)*COS(CSE(J)*PI/180.)/SC 
YSHIP(2)=YS(J) + (SL/4. ) *COS(CSE(J) *PI/180. )/SC + 

* (SB/2. ) *COS( (CSE(J)+90. ) *PI/180. )/SC 
YSHIP(3)=YS(J) - (SL/2. ) *COS(CSE(J) *PI/180. )/SC + 

* (SB/2 . ) *COS ( (CSE (J) +90. ) *PI/180. ) /SC 
YSHIP(4)=YS(J) - (SL/2. ) *COS(CSE(J) *PI/180. )/SC + 

* (SB/2. ) *COS( (CSE(J)-90. ) *PI/180. )/SC 
YSHIP(5)=YS(J) + (SL/4.)*COS(CSE(J)*PI/180.)/SC + 

* (SB/2 . ) *COS ( (CSE ( J) -90 . ) *PI/180 . ) /SC 
C THIS PLOT-88 SUBROUTINE FILLS THE SHIP POLYGON 

C 

, CALL STFILL(l) 

CALL FILL (xship,yship, 5) 

RETURN 

END 

c**************************************************************** 
SUBROUTINE RESCLL(J) 

C THIS SUBROUTINE CALCULATES THE INSTANTANEOUS COORDINATES OF 
C THE RADAR RESOLUTION CELL 

COMMON/MIS S/RG , THETA, AG , PI 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/EQ/ XET(60) , YET(60) 

REAL XCELL(4) , YCELL(4) 

C SET UP THE RESOLUTION CELL OF THE MISSILE 

C 

RMBD = BRGTM(J) + 180. 

AG = RMT(J) * TAN(THETA*PI/180. ) 

ALPHA = 90. + RMBD 

XWl = XET(J) + AG*SIN(ALPHA*PI/180. ) 

YWl = YET(J) + AG * COS (ALPHA* PI/ 180. ) 

XW2 = XET(J) - AG*SIN(ALPHA*PI/180.) 

YW2 = YET(J) - AG*COS(ALPHA*PI/180. ) 
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no ooo o no on 



c 



XCELL(l) 

YCELL(l) 

XCELL(2) 

YCELL(2) 

XCELL(3) 

YCELL(3) 

XCELL(4) 

YCELL(4) 



XWl 

YWl 

XWl 

YWl 

XW2 

YW2 

XW2 

YW2 



RG*SIN 

RG*COS 

RG*SIN 

RG*COS 

RG*SIN 

RG*COS 

RG*SIN 

RG*COS 



(RMBD*PI/180. ) 
(RMBD*PI/180.) 
(RMBD*PI/180.) 
(RMBD*PI/180.) 
(RMBD*PI/180.) 
(RMBD*PI/180. ) 
(RMBD*PI/180.) 
(RMBD*PI/180.) 



DRAW THE RESOLUTION CELL 



CALL STFILL(3) 

CALL FILL (XCELL,YCELL,4) 

RETURN 

END 

*************************************************************** 

SUBROUTINE WIND 

THIS PLOT-88 SUBROUTINE DISPLAYS THE WIND DIRECTION LINE 
COMMON/MISS/RG , THETA, AG , PI 
COMMON/WD/GW , VW 
CHARACTER* 20 LW 
LW ='W' 

XCENT =4.0 
YCENT =3.0 

DRAW WIND DIRECTION LINE 
TWBRG = GW + 180. 

XWDPl = XCENT - 2.8*SIN(TWBRG*PI/180. ) 

YWDPl = YCENT - 2 . 8*COS (TWBRG*PI/180 . ) 

XWDP2 = XCENT - 2 . 0*SIN (TWBRG*PI/180 . ) 

YWDP2 = YCENT - 2 . 0*COS (TWBRG*PI/180 . ) 

CALL PL0T(XWDP1, YWDPl, 3) 

CALL PLOTD(XWDP2,YWDP2,2) 

CALL SYMB0L(XWDP1, YWDPl, 0.20,LW, 0. 0, 1) 

RETURN 

END 

*************************************************************** 

SUBROUTINE CHAFF (J) 

THIS SUBROUTINE CALCULATES THE COORDINATES THE COORDINATES 
OF THE HEXAGON REPRESENTING THE CHAFF 

COMMON/TRACK/ XM (60 ) , YM ( 60) , XS (60) , YS ( 60) , XC ( 60) , YC ( 60 ) 
COMMON/ABCD/BRGM(60) ,R(60) , SC, VS (60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/CFCL/PARTC , CRAD , RGC , AGC 
COMMON/MISS/RG , THETA, AG , PI 
REAL XCF(6) , YCF(6) 

DRAW CHAFF CLOUD 



XCF(1)=XC(J) 

YCF(1)=YC(J) 

XCF(2)=XC(J) 

YCF(2)=YC(J) 

XCF(3)=XC(J) 



+ CRAD 

+ CRAD * 
+ CRAD * 
- CRAD * 



COS(60.*PI/180. ) 
SIN(60.*PI/180. ) 
COS(60.*PI/180. ) 
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YCF(3)=YC(J) + CRAD * SIN(60. *PI/180. ) 

XCF(4)=XC(J) - CRAD 
YCF(4)=YC(J) 

XCF(5)=XC(J) - CRAD * COS (60.*PI/180. ) 

YCF(5)=YC(J) - CRAD * SIN(60. *PI/180 . ) 

XCF(6)=XC(J) + CRAD * COS (60. *PI/180. ) 

YCF(6)=YC(J) - CRAD * SIN(60. *PI/180. ) 

CALL STFILL(2) 

CALL FILL (XCF,YCF,6) 

RETURN 

END 

C**************************************************************** 

SUBROUTINE MISSLE(J) 

C THIS PLOT-88 SUBROUTINE DISPLAYS THE MISSILE BEARING /LINE OF 
FIRE 

COMMON/EQ/ XET(60) ,YET(60) 

COMMON/MISS/RG , THETA , AG, PI 

COMMON/TRACK/ XM(60) ,YM(60) ,XS(60) ,YS(60) ,XC(60) ,YC(60) 
COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/FG/IFLAG1(60) ,IFLAG2(60) 

CHARACTER* 20 LM 
LM = 'M' 

XMl = XET(J) + 2.5 * SIN(BRGTM(J) *PI/180. ) 

YMl = YET(J) + 2.5 * COS ( BRGTM ( J) *PI/180 . ) 

C DRAW THE MISSILE BEARING LINE 

C 

IF (IFLAGl(J) .EQ. 0) THEN 
IF (RMT(J) .GT.2.5) THEN 

CALL PLOT(XET(J) ,YET(J) ,3) 

CALL PL0T(XM1, YM1,2) 

CALL SYMB0L(XM1, YMl, 0.20, LM, 0.0, 1) 

ELSE 

CALL PLOT(XM(J) , YM(J) , 3) 

CALL PLOT(XET(J) , YET(J) ,2) 

CALL SYMBOL (XM(J) , YM ( J) , 0 . 20 , LM, 0 . 0 , 1) 

END IF 
ELSE 

CALL SYMBOL (XM(J) , YM ( J) , 0 . 20 , LM, 0 . 0 , 1 ) 

END IF 
RETURN 
END 

C**************************************************************** 
SUBROUTINE TABLEl(J) 

C THIS PLOT-88 SUBROUTINE DISPLAYS THE VALUES OF IMPORTANT 
VARIABLES 

C BOTH FOR THE INITIAL PICTURE AND FOR THE RUN OF PLAY EVERY 5 
SECS 

COMMON/ TME/T 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/ WD/ GW , VW 
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COMMON/SHCL/AGS , RGS , BRGMS ( 60) , SAA, PARTS , SAAN (60) 
COMMON/RC/RCSS , RCSSP ( 60) , RCSC, RCSCP (60) 

CHARACTER*20 
CS , SPD , WND , MLL, RGE , SCL, TR , AA , SRCS , CRCS , SP, BR, SCR , TM 



TM 


= 




BR 


= 


'BRG' 


SP 




'SHIP' 


CS 


= 


'CSE' 


SPD 


= 


'SPD' 


MLL 




'MISSILE' 


RGE 


= 


'RG' 


WND 


= 


'WIND' 


SCR 


= 


'SCREEN' 


SCL 


= 


'SCALE' 


TR 


= 


'TRUE' 


SRCS 




'SHIP RCS 


CRCS 




'CHAFF RCS 



C CONVERT BACK TO THE UNSCALED RANGE & SPEED 

C 

ZJ = J 

RR = R(J) * SC 

SHSPD = 1.8 * SC * VS(J) 

C DRAW TABLE AT RIGHT OF THE SCREEN 

C 

CALL PLOT (6. 9, 6. 0,3) 

CALL PLOT (6. 9, 0.0, 2) 

call symbol(7.0,5.8,0.18,SP,0.0,4) 
call symbol (7 . 0, 5 . 5 , 0 . 16, CS, 0 . 0, 3) 
call number (7 . 6, 5 .5,0. 16, CSE ( J) , 0 . 0, -1) 
call symbol (7 . 0, 5. 3 , 0 . 16, SPD, 0.0,3) 
call number (7 . 6, 5 . 3 , 0. 16, SHSPD, 0. 0, -1) 
call symbol (7 . 0, 4 . 8 , 0. 16, WND, 0 . 0 , 4 ) 
call symbol (7 . 0, 4 . 5, 0. 16, BR, 0. 0, 3) 
call number (7 . 6, 4 . 5, 0.16, GW, 0. 0,-1) 
call symbol (7. 0,4. 0,0. 16, SPD, 0.0, 3) 
call number (7 . 6, 4 .0,0.16,VW,0.0,-1) 
call symbol (7. 0, 3. 5,0. 16, MLL, 0.0, 7) 
call symbol (7 . 0, 3 . 0, 0. 16, BR, 0. 0, 7) 
call number (7 . 6, 3 . 0, 0 . 16, BRGM ( J) ,0.0,-l) 
call symbol (7. 0, 2 . 5, 0. 16,RGE, 0.0,2) 
call number (7 . 4,2.5,0.16,RR,0.0,-1) 
call symbol (7 . 0, 1. 7 , 0. 16, SCR, 0.0,6) 
call symbol (7 . 2 , 1 . 4 , 0. 16, SCL, 0.0,5) 
call number (7 .3, 1.1, 0.16, SC, 0.0,1) 
call symbol (7 . 1 , 0. 8 , 0. 16,TM, 0. 0, 3) 
call number (7. 4, 0.8,0. 16, ZJ,0. 0,-1) 

C SHIP RCS 

C 

call symbol (4.0,0. 1 , 0 . 16 , SRCS , 0 . 0 , 9) 
call number(6.0,0.1,0.16,RCSSP(J) ,0.0,1) 
C CHAFF RCS 

C 
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call symbol (4.0,0.3,0.16, CRCS ,0.0,9) 
call number (6. 0, 0. 3 ,0. 16,RCSCP(J) ,0.0,1) 

RETURN 

END 

Q* *************************************************************** 

SUBROUTINE SHPCLL(J) 

C THIS SUBROUTINE CALCULATES THE POTION OF THE SHIP RCS IN THE 
CELL 

COMMON/MISS/RG , THETA , AG , PI 

COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/SHCL/AGS , RGS , BRGMS (60), SAA , PARTS , SAAN ( 60 ) 
COMMON/SHP/SL, SB 
C 

IF (RGS.GT. (RG-SL/(2.0*SC) ) ) THEN 
IF (BRGMS (J) .LE. 0.0) THEN 

BR = 90.0 - (BRGMS (J) * 180. /PI) 

ELSE 

BR = 450.0 - (BRGMS (J) * 180. /PI) 

END IF 

CALL ASPECT(J,BR,SAAR) 

S2R = ABS(RGS-RG) 

IF (SAAR. EQ. 90.0) THEN 
SIR = 0.0 
ELSE 

SIR = S2R / ABS(COS(SAAR*PI/180. ) ) 

END IF 
END IF 

IF (AGS.GT. (AG-SL/(2.0*SC) ) ) THEN 
IF (BRGMS (J) .GT.BRGMT(J) ) THEN 
IF (BRGMS (J) .LE. 0.0) THEN 

BA = 180.0 - (BRGMS (J) * 180. /PI) 

ELSE 

BA = 540.0 - (BRGMS (J) * 180. /PI) 

END IF 
ELSE 

IF (BRGMS(J) .LE.0.0) THEN 

BA = 0.0 - (BRGMS(J) * 180. /PI) 

ELSE 

BA = 360.0 - (BRGMS (J) * 180. /PI) 

END IF 
END IF 

CALL ASPECT (J, BA, SAAA) 

S2A = ABS(AGS-AG) 

IF (SAAA. EQ. 90.0) THEN 
SIA = 0.0 
ELSE 

SIA = S2A / ABS(COS(SAAA*PI/180. ) ) 

END IF 
END IF 

IF ( (RGS.GT. (RG-SL/(SC*2.0) ) ) .AND. (AGS.GT. (AG-SL/ (SC*2 . 0) ) ) ) 

THEN 
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IF ( (RGS.LT.RG) .AND. (AGS.LT.AG) ) THEN 
IF (SIA.LT.SIR) THEN 
SI = SIA 
SAA = 90 - SAAA 
ELSE 

SI = SIR 
SAA = SAAR 
END IF 

ELSE IF ( (RGS.GT.RG) .AND. (AGS.GT.AG) ) THEN 
IF (SIA.LT.SIR) THEN 
SI = SIR 
SAA = SAAR 
ELSE 

SI = SIA 

SAA * 90. - SAAA 
END IF 

ELSE IF ( (RGS.GT.RG) .AND. (AGS. LE. AG) ) THEN 
SI = SIR 
SAA = SAAR 
ELSE 

SI = SIA 

SAA = 90.0 - SAAA 
END IF 



ELSEIF( (RGS.GT. (RG-SL/ (SC*2 . 0) ) ) .AND. (AGS.LE. (AG-SL/ (SC*2 . 0) ) ) ) 
*THEN 

SI = SIR 
SAA = SAAR 
ELSE 

SI = SIA 

SAA = 90. - SAAA 
END IF 

SAREA = SL * SB 
SAAN(J) = SAA 
SISC = S1*SC 

IF ( (RGS.GT.RG) .OR. (AGS.GT.AG) ) THEN 

PARTS = (SAREA - .5*SAREA - SB*S1SC) /SAREA 
IF (SISC.GT. (SL/2.0) ) THEN 
PARTS =0.0 
END IF 

ELSE IF (( (RGS.EQ.RG) .AND. (AGS.LE. AG) ) .OR. 

* ( (RGS.LT.RG) .AND. (AGS.EQ.AG) ) ) THEN 
PARTS = 0.5*SAREA 
ELSE 

PARTS = (0.5* SAREA + SB* SISC) /SAREA 
IF (SISC.GT. (SL/2.0) ) THEN 
PARTS =1.0 
END IF 
END IF 
RETURN 
END 
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c**************************************************************** 

SUBROUTINE ASPECT (J, B,SAA) 

C THIS SUBROUTINE CALCULATES THE ANGLE BETWEEN THE SHIP'S HEADING 
C AND THE GATE BOUNDARY IT EXITS. IT IS USED IN SUBROUTINE SHPCLL 
COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

DETERMINE SHIP ASPECT ANGLE (SAA) 



QUADRANT II 



IF ( (90.0. LE.B) .AND. (B.LT. 180.0) ) THEN 

IF ( (0.0. LE. CSE(J) ) .AND. (CSE(J) .LT. 90.0) ) THEN 
IF (B.GE.CSE(J) ) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) + 180.0 - B 
END IF 
END IF 

IF ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) ) THEN 
IF (B.GE.CSE(J) ) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) - B 
END IF 
END IF 

IF ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) ) THEN 
IF ( (CSE(J)-B) .LE.90.0) THEN 
SAA = CSE(J) - B 
ELSE 

SAA = B - (CSE(J) - 180.0) 

END IF 
END IF 

IF ( (270.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 360.0) ) THEN 
IF ((CSE(J)-B) .LE. 180.0) THEN 
SAA = B - (CSE(J) - 180.0) 

ELSE 

SAA = (CSE(J) - 180.0) - B 
END IF 
END IF 
QUADRANT I 



ELSE IF ((180. 0. LE.B) .AND. (B.LT. 270.0) ) THEN 
IF ( (0.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 90.0) ) THEN 
IF (B.GE. (CSE(J)+180) ) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) + 180.0 - B 
END IF 
END IF 

IF ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) ) THEN 
IF ( (B-CSE(J) ) .LE.90. 0) THEN 
SAA = B - CSE(J) 
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ELSE 

SAA = CSE(J) + 180.0 - B 
END IF 
END IF 

IF ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) ) THEN 
IF (CSE(J) .GE.B) THEN 

SAA = CSE(J) - B 
ELSE 

SAA = B - CSE(J) 

END IF 
END IF 

IF ( (270.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 360.0) ) THEN 
IF ( (CSE(J)-B) .LE.90.0) THEN 

SAA * CSE(J) - B 
ELSE 

SAA = B - (CSE(J) - 180.0) 

END IF 
END IF 
QUADRANT IV 



ELSE IF ( (270.0. LE.B) .AND. (B.LT. 360.0) ) THEN 

IF ( (0.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 90.0) ) THEN 
IF ( (B-(CSE(J)+180) ) .LE.90.0) THEN 
SAA = B - (CSE(J)+180. ) 

ELSE 

SAA = CSE(J) - B + 360.0 
END IF 
END IF 

IF ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) ) THEN 
IF (B.GE. (CSE(J)+180.0) ) THEN 
SAA = B - (CSE(J)+180.0) 

ELSE 

SAA = CSE(J) + 180.0 - B 
END IF 
END IF 

IF ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) ) THEN 
IF ( (CSE(J)+180.0-B) .GE.90.0) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) + 180.0 - B 
END IF 
END IF 

IF ((270. O.LE.CSE(J)) .AND. (CSE(J) .LT. 360.0)) THEN 
IF (B.GE.CSE(J) ) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) - B 
END IF 
END IF 

QUADRANT III 



ELSE IF ( (360.0. LE.B) .AND. (B.LT. 450.0) ) THEN 
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IF ( (0.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 90.0) ) THEN 
IF ( (B-360.0) .LE.CSE(J) ) THEN 
SAA = CSE(J) - (B - 360.0) 

ELSE 

SAA = (B - 360.0) - CSE(J) 

END IF 
END IF 

IF ( (90.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 180.0) ) THEN 
IF ( (CSE(J)-(B-360.0) ) .LE.90.0) THEN 
SAA = CSE(J) - (B-360.0) 

ELSE 

SAA = (B - 360.0) - (CSE(J) - 180.0) 

END IF 
END IF 

IF ( (180.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 270.0) ) THEN 
IF ( (CSE(J)-180.0) .GE. (B-360.0) ) THEN 
SAA = (CSE(J) - 180.) - (B - 360.0) 

ELSE 

SAA = (B - 360.) - (CSE(J) - 180.0) 

END IF 
END IF 

IF ( (270.0. LE.CSE(J) ) .AND. (CSE(J) .LT. 360.0) ) THEN 
IF ( (B-CSE(J) ) .LE.90.0) THEN 
SAA = B - CSE(J) 

ELSE 

SAA = CSE(J) - (B - 180.0) 

END IF 
END IF 
ELSE 
END IF 
RETURN 
END 

C *************************************************************** 
SUBROUTINE CHFCLL 

C THIS SUBROUTINE CALCULATES THE CHAFF RCS IN THE CELL 
COMMON/MI SS/RG , THETA , AG , PI 
COMMON/ CFCL/ PARTC , CRAD , RGC , AGC 
COMMON/ RC/RCSS , RCSSP (60), RCSC , RCSCP (60) 

CAREA = PI * CRAD**2.0 
RIR = ABS(RGC - RG) 

RIA = ABS(AGC - AG) 

IF ( (RIR. LT. CRAD) .AND. (RIR. GT. 0.0) ) THEN 
A2R = ACOS( RIR/ CRAD) 

R2R = CRAD*SIN(A2R) 

SECTR = A2R*CRAD**2.0 
TRIR = R1R*R2R 
END IF 

IF ( (RIA. LT. CRAD) .AND. (RIA. GT. 0.0) ) THEN 
A2A = ACOS (RIA/ CRAD) 

R2A = CRAD*SIN(A2A) 

SECTA = A2A*CRAD**2.0 
TRIA = R1A*R2A 
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END IF 

IF ( ( (RIR.LT.CRAD) .AND. (RlR.GT.0.0) ) .AND. 

*( (RIA.LT.CRAD) .AND. (RlA.GT.0.0) ) ) THEN 
A2A = ACOS (RIA/CRAD) 

A2R = ACOS (RIR/CRAD) 

R2A = CRAD*SIN(A2A) 

R2R = CRAD*SIN(A2R) 

SECTA = A2A*CRAD**2 . 0 
SECTR = A2R*CRAD**2 . 0 
TRIA = R1A*R2A 
TRIR = R1R*R2R 
END IF 

IF ( (RlR.EQ.0.0) .OR. (RlA.EQ.0.0) ) THEN 
PARTC =0.5 

IF ( (RlR.EQ.0.0) .AND. (RlA.EQ.0.0) ) THEN 
PARTC =0.25 
END IF 

ELSE IF ( (RIR.LT.CRAD) .AND. (RIA.GT.CRAD) ) THEN 
IF (RGC.LT.RG) THEN 

PARTC = (CAREA-SECTR+TRIR)/CAREA 
ELSE 

PARTC = (SECTR-TRIR)/CAREA 
END IF 

ELSE IF ( (RIR.GT.CRAD) .AND. (RIA.LT.CRAD) ) THEN 
IF (AGC.LT.AG) THEN 

PARTC = (CAREA-SECTA+TRIA)/CAREA 
ELSE 

PARTC = (SECTA-TRIA)/CAREA 
END IF 

ELSE IF (( (RIR.LT.CRAD) .AND. (RIR.GT. 0.0) ) .AND. 

*( (RIA.LT.CRAD) .AND. (RlA.GT.0.0) ) ) THEN 
IF ( (RGC.LT.RG) .AND. (AGC.LT.AG) ) THEN 
IF (RIR.LT.RIA) THEN 

PARTC = (CAREA-SECTR+TRIR)/CAREA 
ELSE 

PARTC = (CAREA-SECTA+TRIA)/CAREA 
END IF 

ELSE IF ( (RGC.LT.RG) .AND. (AGC.GT. AG) ) THEN 
PARTC = (SECTA-TRIA)/CAREA 
ELSE IF ( (RGC.GT.RG) .AND. (AGC.LT.AG) ) THEN 
PARTC = (SECTR-TRIR)/CAREA 
ELSE 

IF (RIR.LT.RIA) THEN 

PARTC = 0.5*(SECTA-TRIA)/CAREA 
ELSE 

PARTC = 0.5*(SECTR-TRIR)/CAREA 
END IF 
END IF 
ELSE 

PARTC =0.0 
END IF 

IF ( (RGC.GT. (RG+CRAD) ) .OR. (AGC.GT. (AG+CRAD) ) ) THEN 



111 



PARTC =0.0 
END IF 
RETURN 
END 

C* *************************************************************** 

SUBROUTINE SCREN3 (IFIAG3) 

C THIS SUBROUTINE ENABLES ALTERATION OF THE INPUT VARIABLES 
COMMON/ABCD/BRGM(60) ,R(60) ,SC,VS(60) ,CSE(60) ,BRGMT(60) , 
*BRGTM(60) ,RMT(60) 

COMMON/WD/GW , VW 
COMMON/MISS/RG , THETA, AG , PI 
COMMON/RC/RCSS , RCSSP ( 60 ) , RCSC , RCSCP ( 60 ) 

INTEGER I 



IFLAG3=0 
100 WRITE (*,*) 
WRITE(*,*) 
WRITE (*, *) 
WRITE(*, *) 
WRITE(*,*) 
WRITE(*, *) 
WRITE(*,*) 
WRITE(*,*) 
WRITE(*, *) 
WRITE (*, *) 
WRITE (*, *) 
WRITE (*,*) 
WRITE(*,*) 
WRITE (*,*) 
READ (*,*) 
IF (I.EQ.l) 



( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

( 9 ) 

( 10 ) 
( 11 ) 

'( 12 ) 

' (13) 

' INPUT 
I 

THEN 



MISSILE BEARING' 

MISSILE RANGE' 

SHIP"S HEAD' 

SCREEN SCALE' 

MISSILE RADAR BEAMWIDTH' 

RANGE GATE' 

MAXIMUM RADAR CROSS SECTION OF 
MINIMUM RADAR CROSS SECTION OF 
RADAR CROSS SECTION OF CHAFF' 
TRUE WIND DIRECTION' 

TRUE WIND SPEED' 

ALL PARAMETERS' 

RETURN TO MAIN PROGRAM' 
SELECTION' 



SHIP' 

SHIP' 



WRITE(*, *) 

READ (*,*) 

IFLAG3 = 2 
GO TO 100 
ENDIF 

IF (I.EQ.2) THEN 
WRITE (*,*) 'INPUT 
READ (*,*) RANGE 
IFLAG3 = 2 
GO TO 100 
ENDIF 

IF (I.EQ.3) THEN 
WRITE (*,*) 'INPUT 
READ (*,*) CSE(l) 
IFLAG3 = 2 
GO TO 100 
ENDIF 

IF (I.EQ.4) THEN 
WRITE (*,*) 'INPUT 
READ (*,*) SC 
IFLAG3 = 2 
GO TO 100 



'INPUT MISSILE 
BRGM(l) 



BEARING' 



INITIAL MISSILE RANGE' 



SHIP"S HEAD' 



SCREEN SCALE' 
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ENDIF 

IF (I.EQ.5) 
WRITE(*,*) 
READ (*,*) 
IFLAG3 = 2 
GO TO 100 
ENDIF 

IF (I.EQ.6) 
WRITE (*,*) 
READ (*,*) 
RG=RG/SC 
IFLAG3 = 2 
GO TO 100 
ENDIF 

IF (I.EQ.7) 
WRITE (*,*) 
READ (*,*) 
I FLAGS = 2 
GO TO 100 
ENDIF 

IF (I.EQ.8) 
WRITE (*, *) 
READ (*,*) 
I FLAGS = 2 
GO TO 100 
ENDIF 

IF (I.EQ.9) 
WRITE (*, *) 
READ (*,*) 
IFLAGS = 2 
GO TO 100 
ENDIF 

IF (I.EQ.IO) 
WRITE (*,*) 
READ (*,*) 
IFLAGS = 2 
GO TO 100 
ENDIF 

IF (I.EQ.ll) 
WRITE(*,*) 
READ (*,*) 
IFLAGS = 2 
GO TO 100 
ENDIF 

IF (I.EQ.12) 
IF ((I.LT.l) 
WRITE {*,*) 
GO TO 100 
ENDIF 
RETURN 
END 



THEN 

'INPUT MISSILE RADAR BEAMWIDTH' 
THETA 



THEN 

'INPUT RANGE GATE' 
RGATE 



THEN 

'INPUT MAXIMUM RADAR CROSS SECTION OF SHIP' 
RCSSMAX 



THEN 

'INPUT MINIMUM RADAR CROSS SECTION OF SHIP' 
RCSSMIN 



THEN 

'INPUT RADAR CROSS SECTIONOF CHAFFf' 
RCSSC 



THEN 

'INPUT TRUE WIND DIRECTION' 
GW 



THEN 

'INPUT TRUE WIND SPEED' 
VW 



IFLAGS=1 

OR. (I.GT. IS) ) THEN 
'INPUT AGAIN' 



IIS 
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